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HIGH  TEMPERATURE  POLYMER  PROGRAM 


CONDENSATION  REACTIONS  OF  HEXAPHENYL- 
DICHLOROPHOSPHO NITRILE  TETRAMER 
WITH  DIOLS  AND  DIAMINES 

by 

D.  L.  Herring,  C.  M.  Douglas,  and  A.  J.  Bilbo 
Chemistry  Division 


INTRODUCTION 

The  synthesis  of  hexaphenyldichlorophosphonitrile  tetramer,  compound 
(1 ),  by  the  reaction 


Cl 


(C6H5)2PC1  +  C6H5PC12  +  NaN3‘ 


■n 


where  I  I  represents 


He, 


(I) 


f6H5  ^6H5 

—  P — N=P — C,  H_ 
i  o  5 


N  N 


C ,  H_  —  P=N —  P — 

6  5  I  l 

C6H5  C6H5 

has  been  described  in  an  earlier  report.^  In  addition  to  the  diethoxy* 
and  diphenoxy-derivatives  that  were  also  reported,  further  substitution 
reactions  have  been  carried  out  and  the  dipiperidino-,  diamino-,  and 
bis(dimethylamino)-hexaphenylphosphonitrile  tetramer s— compounds 
(II),  (III),  and  (IV),  respectively— have  been  prepared. 


*D.  L.  Herring  and  C.  M.  Douglas,  NAVWEPS  Report  8226,  Quarterly 
Report;  Foundational  Research  Projects,  October- December  1964,  p.  5. 


DISCUSSION 


Upon  pyrolysis,  compound  (III)  condenses,  with  the  elimination  of 
ammonia,  to  form  a  low  molecular  weight  species  which  is  probably 
dimeric: 


H 

N 


H 


Heating  compound  (I)  at  265 *C  in  a  sealed  tube  results  in  the  forma¬ 
tion  of  a  glassy  material  (m.p.  70-80*C).  By  fractional  crystallization, 
the  crude  product  can  be  separated  into  the  isomers  of  hexaphenyl- 
dichlorophosphonitrile  tetramer.  Compound  (I),  m.p.  303-304*C,  is 
recovered  in  30%  yield,  and  the  remainder,  compound  (V),  melts  at 


185-187*C. 


The  thermal  rearrangement  of  compound  (I)  was  accompanied  by  the 
formation  of  a  very  small  quantity  of  benzene,  that  is,  1  mmole  of 
benzene  per  278  mmoles  of  compound  (I).  It  can  be  postulated  that  the 
benzene  was  formed  by  condensation  of  2  molecules  of  (I),  thus: 


Cl  _ 

2  1 


+  C6H6 


(VI) 


Some  support  for  the  above  explanation  is  given  by  consideration  of  the 
properties  of  the  polymeric  products  obtained  from  the  condensation  of 
the  isomer  mix  and  biphenol.  The  melting  points  were  above  300*C  and 
benzene  insolubility  was  observed.  These  are  properties  of  cross- 
linked  material  and  would  be  expected  if  compound  (VI)  was  present  in 
the  mixture  of  isomers. 
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It  is  believed  that  the  higher-melting  isomer  has  the  trans  configu¬ 
ration  and  that  compound  (V)  has  the  cis  configuration.  Compound  (V), 
when  treated  with  LiN3  in  refluxing  acetonitrile,  yields  a  mixture  of 
cis-  and  trans-diazido  tetramers,  as  docs  compound  (I).  A  careful  study 
of  the  dipole  moments  of  low-  and  high-melting  diazido  tetramers  has 
proven  them  to  be  the  cis  and  trans  forms,  respectively. 2 

Compound  (I)  represents  a  useful  starting  material  for  the  preparation 
of  polymers  in  which  the  phosphonitrile  ring  system  can  be  linked  in  a 
systematic  linear  fashion  by  means  of  other  difunctional  compounds  such 
as  diols. 


+  HOROH 
Cl 


For  the  polymerization  studies,  compound  (I)  was  recrystallized  twice 
from  toluene,  dried  at  150°C  for  24  hours,  and  stored  in  evacuated  tubes. 
Biphenol,  which  was  supplied  by  Eastman  Organic,  was  sublimed 
twice.  It  was  then  finally  sublimed,  and  only  the  middle  80%  was 
retained. 

The  ratio  of  (I)  to  biphenol  was  varied  from  1:1.2  to  1:0.95  to  deter¬ 
mine  the  effect  of  any  inert  impurities  that  might  be  present  in  the 
starting  materials.  This  imbalance  in  the  stoichiometry  had  no  appre¬ 
ciable  influence  either  on  the  range  of  molecular  weights  (10,000  to 
12,000)  or  in  the  softening  points  (185-220*C).  In  all  cases,  the  evolu¬ 
tion  of  hydrogen  chloride  gas  was  over  95%  of  theoretical.  However, 
elemental  analysis  of  most  of  the  polymeric  resins  showed  no  chlorine. 

A  temperature  of  290-220 *C  was  required  for  the  initiation  of  the 
melt  reaction,  which  proceeded  very  slowly  at  this  temperature,  but 
continuation  was  aided  by  the  removal  of  the  hydrogen  chloride  generated. 
Although  it  was  convenient  to  carry  out  the  polymerisations  at  reduced 
pressure,  it  was  noted  that  at  220 *C  biphenol  lublimed  out  of  the  reac¬ 
tion  zone.  When  the  bath  temperature  was  raised,  additional  diol  was 
lost,  making  a  total  loss  up  to  2%.  Carrying  out  the  polymerisation 
under  nitrogen  or  in  sealed  tubes  that  were  relieved  of  HC1  pressure 
periodically  did  not  improve  the  results. 

It  was  anticipated  that  use  of  the  lower-melting  isomer,  compound  (V), 
with  biphenol  would  give  a  mixture  which  would  react  at  such  a  low  tem¬ 
perature  that  loss  of  the  diol  would  not  occur.  However,  even  though 


2 


C.  M.  Sharts  and  O.  Gentry,  unpublished  results. 
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the  reaction  mixture  was  held  at  the  melting  point  of  (V),  no  gas  was 
evolved  and  prolonged  heating  caused  the  sublimation  of  the  diol. 

As  pointed  out  earlier,  the  unseparated  mixture  of  the  phosphonitrilic 
isomers  melts  at  70-80*C.  The  polymerization  of  this  mixture  with 
biphenol  yielded  a  resinous  material  which  softened  slightly  at  320 *C  but 
did  not  melt.  It  was  insoluble  in  most  organic  solvents  and  was  swollen 
by  benzene.  No  molecular  weight  measurements  could  be  made. 

Methods  for  circumventing  the  loss  of  the  diol  from  the  reaction  zone 
have  also  included  the  use  of  the  disodium  salt  of  biphenol.  Results 
from  one  preliminary  experiment  have  been  encouraging  in  that  a  resin 
which  did  not  melt  at  320 *C  was  isolated.  Again,  however,  insolubility 
of  this  product  prevented  the  molecular  weight  determination. 

Solution  polymerizations  were  carried  out  in  refluxing  anisole, 
diphenyl  ether,  petrolatum,  and  acetonitrile  containing  quinoline.  Hydro¬ 
gen  chloride  was  slowly  generated  over  a  six-week  period  in  the  anisole 
reaction.  After  removal  of  the  solvent,  further  heating  of  the  product 
in  vacuum  not  only  yielded  additional  HCl,  but  resulted  in  the  sublima¬ 
tion  of  unreacted  diol.  The  polymer  was  a  low  molecular  weight  species 
(5200)  containing  no  chlorine  as  end  groups.  The  reaction  in  the  higher 
boiling  solvent,  diphenyl  ether,  proceeded  more  rapidly,  but  did  not  near 
completion  after  a  week  of  refluxing.  About  75%  of  the  theoretical  amount 
of  HCl  was  obtained.  A  mixture  of  (1)  and  biphenol  suspended  in  boiling 
petrolatum  yielded  only  a  small  fraction  of  the  expected  HCl.  An  equi¬ 
molar  mixture  of  (1)  and  biphenol  did  not  react  in  refluxing  acetonitrile 
with  or  without  the  presence  of  an  acid  acceptor  such  as  quinoline. 

The  compound  (l)-biphenol  polymer  can  be  fractionated,  by  use  of  a 
benzene-heptane  mixture,  into  a  material  having  a  molecular  weight  of 
18,000.  The  melting  range  of  this  fraction  was  not  substantially  higher 
than  the  original  polymer  with  a  molecular  weight  of  M 2,000.  Separa¬ 
tion  of  the  polymer  species  gave,  as  the  most  soluble  portion,  a  crystal¬ 
line  material  whtfe  molecular  weight  and  elemental  analysis  correspond 
very  closely  with  that  calculated  for  the  dimer: 
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Compound  {I)  couples  with  diamines  according  to  the  following  equation: 


Cl 


n 

Cl 


R*  R’ 

I  I 

+  H— N— R— N— H- 


•la 


N — R — N4~ 


+  2nHCl 


n 


where  R'  =  H  or  C.H 

o  5 

R  =  biphenylene 


In  this  condensation,  it  is  possible  that  some  of  the  hydrogen  chloride 
reacts  with  the  amine  to  form  the  hydrochloride.  Since  no  sublimation 
was  noted,  any  amine  hydrochloride  produced  initially  was  probably 
decomposed  at  the  elevated  reaction  temperature.  It  is  interesting  that 
although  the  gas  evolution  indicates  that  the  reaction  was  only  83%  com* 
plete,  the  resinous  product  contains  no  chlorine  and  the  elemental  analy- 
sis  is  very  close  to  that  calculated  for  the  polymer  shown  in  the  brackets 
above  (where  n  is  large). 

When  R'  was  phenyl,  the  coupling  proceeded  less  satisfactorily. 

A  temperature  100*C  higher  was  required  to  initiate  the  reaction,  and 
although  somewhat  more  HC1  was  generated,  the  product  was  a  low- 
melting.  low  molecular  weight  material  containing  0.5%  residual  chlorine. 


Compound  (I)  reacted  with  trialkyl  phosphates  in  the  following  manner: 


♦  (RO)3PO- 


IQ  ' 


OR 


OP — 04- 

II 

o 


+  2nRCl 


In  the  case  of  ethyl  phosphate,  evolution  of  ethyl  chloride  indicated 
that  the  reaction  was  95%  complete.  The  product  was  a  low-melting 
resin  which  was  insoluble  in  benzene. 


In  view  of  the  relatively  high  temperatures  that  must  be  employed  to 
effect  the  tetramer  dichloririe-diol  condensations,  it  appeared  that  it 
would  be  advantageous  to  prepare  the  tetramer  dibromide,  which  might 
be  reactive  at  a  lower  temperature.  Replacement  of  chloride  by  aside 
groups  takes  place  readily  when  the  tetramer  dichloride  and  excess 
lithium  azide  are  allowed  to  react  in  boiling  acetonitrile. ^  It  therefore 


>NAVWEPS  Report  8226,  p.  23. 
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seemed  that  lithium  bromide,  which  is  very  soluble  in  acetonitrile,  might 
engage  in  the  same  type  of  repl«,_ement  reaction.  Despite  prolonged 
contact  between  lithium  bromide  and  tetramer  dichloride  in  refluxing 
acetonitrile,  no  reaction  occurred  and  the  tetramer  dichloride  was 
recovered  unchanged. 


EXPERIMENTAL 


Preparation  of  Hexaphenyldipiperidino- 
phosphonitrile  Tetramer  (II) 

Hexaphenyldichlorophosphonitrile  tetramer  (3.0  g)  was  dissolved  in 
100  ml  oi  refluxing  chlorobenzene,  and  an  excess  (10  cc)  of  piperidine 
was  added.  Upon  cooling,  a  precipitate  of  piperidine  hydrochloride  was 
formed  and  was  removed  by  filtration.  The  chlorobenzene  solution  was 
concentrated  to  50  ml,  yielding  the  crude  product.  The  dipiperidino- 
derivative  was  recrystallized  from  benzene  to  give  2.0  g  of  white  crys¬ 
talline  material  melting  at  290-292*C.  The  infrared  spectrum  is  shown 
in  Figure  1. 


Anal.  Calc,  for 

Found: 


C.  68.1;  H,  6.2;  P.  15.3;  N,  10.4; 
mol.  wt.  810. 

C,  68.1;  H,  6.4;  P,  15.1;  N,  10.5; 
mol.  wt.  832. 


Preparation  of  Hexaphenyldiamino- 
phosphonitrile  Tetramer  (III) 

A  slurry  containing  20.0  g  of  (I)  in  500  ml  of  chloroform  was  treated 
with  gaseous  ammonia  for  6  hours.  The  reaction  mixture  was  then 
filtered  to  remove  ammonium  chloride.  Evaporation  of  the  chloroform 
yielded  the  crude  product,  which  recrystallized  from  heptane  as  long 
white  needles,  m.p.  205-207*C.  Its  infrared  spectrum  is  shown  in  Fig 
ure  2. 


Anal.  Calc,  for  C-.H.  .P.N.; 
■ . .  36  34  4  6 

Found: 


C,  64.1;  H,  5.0;  P,  18.4;  N,  12.5; 
mol.  wt.  674. 

C,  64.0;  H.  4.9;  P.  18.1;  N,  12.6; 
mol.  wt.  668. 


Pyrolysis  of  Hexaphenyldiamlnophosphonitrile  Tetramer 

A  sample  of  5.0  g  (0.007  mole)  of  compound  (Ill)  was  heated  in  vacuum 
at  25-360*C  for  12  hours.  During  the  heating  cycle  0.006  mole  of  ammonia 
was  generated,  along  with  a  small  amount  of  benzene.  Both  volatile 


I 


FIGURE  2.  Infrared  Spectrum  of  Hexaphenyldiaminophosphonitriie  Tetramer. 

Compound  (HI) 


products  were  identified  by  infrared  analysis.  The  residue  was  a  tan 
resin  softening  over  the  range  160-185 *C. 

Anal.  Calc,  for  C,,HL .  P.N_:  C,  65.8;  H,  4.7;  P,  18.9;  N,  10.7; 

.  JO  jl  4  o  , 

mol.  wt.  per  unit  657. 

Found;  C,  65.6;  H,  4.8;  P,  18.9;  N,  10.8; 
mol.  wt.  1320. 


Preparation  of  Hexaphenylbis(dimethylamino)- 
phosphonitrile  Tetramer  (IV) 


Gaseous  dimethylamine  was  combined  with  a  slurry  of  250  ml  of 
chloroform  and  14.0  g  of  (I),  The  reaction  mixture  was  stirred  vigor¬ 
ously  during  the  addition  of  (CHj^NH.  /iter  1  hour,  the  reaction  was 
complete,  as  indicated  by  the  disappearance  of  (C£,H5)^Cl2p4N4.  The 
chloroform  solution  was  completely  clear.  The  solvent  was  then  evapo¬ 
rated  and  the  solids  were  extracted  with  benzene  to  separate  the  dimethyl- 
amine  hydrochloride.  The  benzene  solution  was  concentrated,  and  when 
it  had  cooled  the  hexaphenylbis(dimethylamino)phosphonitrile  was  re¬ 
covered  as  white  needles,  m.p.  189-19Q°C.  The  infrared  spectrum 
is  shown  in  Figure  3. 


Anal.  Calc,  for  C..H..P.N.: 
-  40  42  4  6 

Found: 


C,  65.8;  H,  5.8;  P,  17.0;  N,  11.5; 
mol,  wt.  730. 

C,  66.0;  H,  6.0:  P,  16.7;  N,  11.6; 
mol.  wt.  723. 


Thermal  Arrangement  of  Hexaphenyldichloro- 
phosphonitrile  Tetramer  (I) 

■} 

A  30-g  sample  of  (I),  which  had  been  recrystallized  twice  from 
toluene,  was  sealed  in  an  evacuated  glass  ampoule  and  heated  in  an  oven 
at  265*C  for  48  hours.  After  cooling,  the  ampoule  was  opened  into  a 
vacuum  line  and  the  volatile  material  was  condensed  into  a  measuring 
system*  Benzene,  3.4  cc  STP  (0.151  mmole),  was  collected  and  identi¬ 
fied  by  its  infrared  spectrum.  The  non-volatile  material,  m.p.  70-80*C, 
was  then  extracted  with  benzene,  and  the  insoluble  portion  (about  25%) 
was  identified  by  its  infrared  spectrum  and  melting  point  (298- 300  *C) 
as  the  starting  material.  The  benzene  was  evaporated  and  the  residue 
was  twice  recrystallized  from  heptane.  It  yielded  an  isomer,  (V),  m.p. 
185-187*C,  with  the  infrared  spectrum  shown  in  Figure  4. 


FIGURE  3.  Infrared  Spectrum  of  Hexaphenylbis(dimethylamino)phosphonitrile  Tetramer, 

Compound  (IV) 
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Reaction  of  Compound  (V)  with  Lithium  Azide 

Two  separate  experiments  were  carried  out.  In  the  first  experiment, 
2.0  g  of  (V)  was  permitted  to  react  with  3.0  g  of  lithium  azide  in  50  ml 
of  refluxing  acetonitrile.  After  24  hours,  solids  were  filtered  off  and 
solvents  removed  until  only  20  ml  of  solution  remained.  A  first  crop 
of  0.18  g  of  crystals  was  obtained,  m.p.  178-212eC»  which  was  identified 
by  infrared  spectroscopy  as  trans -2, 6-diazidocyclophospho nitrile. 

A  second  crop  of  the  same  diazide,  amounting  to  0.10  g,  was  obtained, 
making  a  total  of  0.28  g  (about  15%  of  theory)  of  diazide.  No  additional 
solid  could  be  obtained  from  the  mother  liquor  of  the  crystallizations. 
The  viscous  glass  had  an  infrared  spectrum  which  suggested  the  glass 
was  a  mixture  of  cis-  and  trans-2,6-diazidocyclophosphonitrile.  In 
addition,  the  bands  of  a  hydroxyl  group  were  observed. 

When  the  experiment  was  repeated,  the  same  results  were  obtained. 
An  analysis  of  the  reaction  product  that  could  not  be  crystallized  showed 
it  to  be  free  of  chlorine. 

Anal.  Calc,  for  C0,H,nP„N  •  C,  59.51;  H,  4.16;  N,  1 9.2.8;  P,  17.05. 

-  36  30  4  10 

Found:  C,  60.36;  H,  4.49;  N,  19.35;  P,  15.89. 

Polymerization  of  Hexaphenyldichlorophos- 
phonitrile  Tetramer  and  Biphenol 

MELT  REACTIONS 

A  mixture  of  7.1363  g  (0.01  mole)  of  compound  (I)  and  2.0507  g  (0.011 
mole)  of  biphenol  was  heated  in  an  evacuated  ampoule.  Evolution  of 
hydrogen  chloride  began  at  220-225 °C,  although  the  mixture  was  not 
fused  at  this  temperature.  Over  the  course  of  10  hours,  the  bath  tem¬ 
perature  was  gradually  increased  to  400 #C.  The  small  amount  of  mate¬ 
rial  which  had  sublimed  into  the  neck  of  the  ampoule  during  the  heating 
cycle  was  recovered,  weighed  (0.1145  g),  and  identified  as  biphenol  by 
its  infrared  spectrum.  The  total  amount  of  HC1  generated  was  96.7%  of 
the  calculated  amount.  The  polymeric  product  was  a  transparent  glass 
melting  over  the  range  1 85-220 *C  and  having  a  molecular  weight  of 
11,360. 

A  summary  of  information  on  the  melt  polymerizations  involving 
tetramer  dichloride  and  biphenol  is  shown  in  Table  1. 


This  reaction  was  conducted  by  C.  M.  Sharts. 


TABLE  1.  Polymerization  of  Tetramer  Dichloride 

and  Biphenol 


Molar  Ratio  of 
Tet.  Dichloride 
and  Biphenol 

Species  of 

Tetramer  Dichloride 
Used 

Percent 
of  HC1 
Obtained 

Mol.  Wt. 
of 

Polymer 

Melting 

Range, 

•c 

1.0  : 

:  0.95 

Compound  (I) 

71 

6,420 

165-175 

1.0  :  1.0 

Compound  (I) 

92.7 

10,000 

170-200 

1.0  ; 

:  1.1 

Compound  (I) 

96.7 

11,360 

185-220 

1.0  :  1.2 

Compound  (I) 

98.4 

12,000 

165-185 

1.0  :  1.0 

Low-melting  isomer, 
(V),  185-187#C 

97.0 

14,500 

180-240 

1.0  :  1.0 

Mixture  of  isomers 

98.5 

Insoluble 

>300 

1.0  : 

:  0.95 

Mixture  of  isomers 

96.2 

Insoluble 

>300 

Reaction  of  Compound  (1)  with  Biphenol  in  a  Sealed  Tube.  An  inti¬ 
mate  mixture  of  0.579  g  (0.813  mmole)  of  compound  (I)  and  0.2088  g 
(1.12  mmole)  of  biphenol  was  loaded  into  a  50-cc  reaction  tube.  After 
being  evacuated,  the  tube  was  sealed  off  and  heated  in  an  oven  at  265  *C 
for  8  hours.  The  tube  neck  was  broken  by  a  tube  breaker  attached  to 
the  vacuum  line  and  HC1  gas  was  condensed  into  a  measuring  system. 

The  amount  of  HC1  calculated  for  0.813  mmole  of  (I)  was  36.4  cc  STP; 
found.  15.7  cc  STP,  or  43%.  The  tube  neck  was  cleaned  and  sealed  off 
under  vacuum  a  second  time  and  heated  at  265 *C  for  48  hours.  Only 
1.5  cc  STP  of  HC1  was  found.  The  experiment  was  abandoned  after  sev¬ 
eral  more  cycles  of  heating  and  venting  of  gas  yielded  negligible  amounts 
of  hydrogen  chloride. 

Preparation  of  the  Disodium  Salt  of  Biphenol,  NaO~-<~1~,V-~<  ONa. 

Twenty-five  grams  (0.13  mole)  of  biphenol  was  dissolved  in  200  ml  of 
3.9  molar  sodium  hydroxide.  The  addition  of  800  ml  of  acetone  to  the 
solution  resulted  in  the  precipitation  of  the  disodium  biphenolate.  The 
product  was  recrystalliced  twice  from  400-ml  portions  of  ethanol,  then 
dried  in  vacuum  at  180*C. 

Anal.  Calc,  for  C12H8°2N*2:  C,  62.6;  H,  3.5;  O,  13.9;  Na.  20.0. 

Found:  C,  62.2;  H,  3.7;  O,  13,6;  Na,  20.2. 

Reaction  of  Compound  (I)  and  Disodium  Biphenolate.  In  the  drybox, 
a  mixture  of  5.0217  g  (0.007  mole)  of  compound  (I)  and  1.6219  g  (0.007 


mole)  of  sodium  biphenolate  was  loaded  into  a  150-ml,  round-bottomed 
flask  containing  stainless  steel  balls.  The  flask  was  then  attached  to  a 
Rinco  Evaporator,  slowly  rotated  under  vacuum,  and  heated  to  310#C 
over  a  period  of  24  hours.  After  the  dark  brown  reaction  mixture  was 
extracted  with  boiling  benzene,  the  benzene-soluble  portion  was  dried 
in  vacuum;  its  melting  range  was  150-165aC. 


Anal.  Calc,  for  C^gH^gP^^O^: 

Found: 


C,  69.7;  H,  4.6;  P,  15.0;  N, 
O.  3.9. 

C,  69.2;  H,  4.7;  P,  15.0;  N, 
O,  3.5;  Na,  0.77;  Cl,  1.2. 


6.8; 

6.4; 


The  benzene-insoiuble  fraction  was  extracted  with  boiling  water,  then 
washed  with  ethanol  and  ether,  and  dried.  This  material  softened 
slightly  at  300 *C  but  did  not  melt. 

SOLUTION  POLYMERIZATION 


Reaction  of  Compound  (I)  and  Biohenol  in  Anisole.  Compound  (I) 
(14.26  g,  0.02  mole)  and  3.72  g  (0.02  mole)  of  biphenol  were  dissolved 
in  500  ml  of  anisole,  b.p.  152-154*C.  The  solution  was  refluxed  under 
nitrogen.  After  an  induction  period  of  about  3  hours,  the  evolution  of 
hydrogen  chloride  began.  Heating  was  continued  for  6  weeks  until  no 
additional  HC1  was  detected.  The  anisole  was  then  evaporated  and  the 
solid  product  was  heated  in  vacuum  at  300 *C  for  about  8  hours.  During 
this  heating  cycle,  a  small  amount  of  hydrogen  chloride  w?s  generated 
and  a  few  milligrams  of  unreacted  diol  was  sublimed;  both  of  these  com¬ 
pounds  were  identified  by  infrared  spectra.  The  polymeric  product  was 
a  brittle  brown  glass  softening  over  the  range  140-180*C. 

Anal.  Calc,  for  C  *aH  P  N.O.:  C,  69.7;  H,  4.6;  P,  15.0;  N.  6.8; 

O,  3.9;  mol.  wt.  per  unit,  826. 

Found:  C,  70.0;  H,  4.9;  P,  14.8;  N,  6.6; 

O,  4.2;  Cl,  0;  mol.  wt.  5200. 

Reaction  of  Compound  (I)  and  Biphenol  in  Diphenyl  Ether.  A  mixture 
of  7.13  g  (0.01  mole)  of  compound  (I)  and  1.90  g  (0.0102  mole)  of  b; phenol, 
together  with  50  ml  of  diphenyl  ether  and  a  magnetic  stirrer  bar,  was 
placed  in  a  cylindrical  tube  (40  X  250  mm)  which  had  been  carefully  dried 
and  filled  with  dry  Nj  gas.  The  tube  was  equipped  with  a  removable 
head  which  carried  a  water  condenser,  stopcock,  and  joint  for  attach¬ 
ment  to  the  vacuum  line.  The  apparatus  was  connected  to  a  vacuum  line 
which  had  incorporated  into  it  a  manostat  and  a  manometer  that  was  open 
to  the  atmosphere.  The  system  and  reaction  apparatus  were  filled  with 
dry  N£  at  atmospheric  pressure  and  the  manostat  was  set  to  maintain 
this  pressure. 


An  oil  bath  heated  with  a  hot  plate-magnetic  stirrer  combination  was 
used  to  bring  the  stirred  reaction  mixture  to  the  reflux  temperature  of 
diphenyl  ether,  258-259°C.  As  pressure  in  excess  of  1  atmosphere 
developed  because  of  increasing  diphenyl  ether  vapor  pressure  and  evolu¬ 
tion  of  HC1  gas,  the  manostat  released  gas,  which  passed  through  sev¬ 
eral  liquid- nitrogen- cooled  traps  and  thence  through  the  pumps.  After 
the  material  had  refluxed  for  8  hours,  the  apparatus  was  cooled  and 
evacuated.  Volatile  components  were  caught  in  the  above-mentioned 
traps.  After  fractionation  of  the  volatiles,  252.8  cc  STP  (56.4%)  of  'HCl 
gas  was  collected.  The  amount  of  HC1  calculated  for  0.01  mole  of  (I) 
was  444  cc  STP.  The  above  treatment  was  repeated  and  continued  for 
6  days  but  the  total  HC1  evolved  amounted  to  only  347  cc  STP,  or  77.6% 
of  theory. 

The  reaction  mixture  was  not  worked  up  because  of  the  low  yield  of 
HC1,  which  indicates  a  poor  condensation  reaction. 

Attempted  Reaction  of  (I)  and  Biphenol  in  Petrolatum.  An  equimolar 
mixture  of  compound  (I)  and  biphenol  suspended  in  distilled  petrolatum 
was  heated  and  stirred  at  275 *C  in  the  apparatus  described  above.  Only 
a  small  fraction  of  the  theoretical  yield  of  HC1  gas  was  obtained. 

Attempted  Reaction  of  Compound  (1)  and  Biphenol  in  Acetonitrile 
with  Quinoline  as  Acid  Acceptor.  In  the  drybox,  2  g  of  an  equimolar 
mixture  of  compound  (I)  and  biphenol  was  added  to  a  flask  containing 
250  cc  of  dry  acetonitrile  and  1.3  g  (0.0 J  mole)  of  freshly  distilled 
quinoline.  After  the  mixture  was  refluxed  for  48  hours  and  filtered, 

1.53  g  (95.6%)  of  compound  (I)  was  recovered  unchanged.  Identification 
was  made  by  infrared  spectrum. 

Fractionation  of  Compound  (l)-Biphenol  Polymer 

The  polymer  sample  used  for  this  preliminary  fractionation  was  pre¬ 
pared  by  the  melt  reaction  described  earlier  (page  12).  The  molecular 
weight,  determined  osmometrically,  was  11,360.  Bensene  was  used 
as  the  solvent  and  n-heptane  as  the  non- solvent.  A  summary  of  the  data 
obtained  from  the  fractionation  of  a  3.9474-g  sample  is  shown  in  Table  2. 

The  fourth  fraction  was  a  crystalline  material  with  an  elemental 
analysis  corresponding  closely  to  that  calculated  for  C48H38^4^402* 

The  infrared  specjrum  (Figure  5)  was  characterised  by  absorption  bands 
much  sharper  than  those  encountered  in  the  tetramer  dichloride-biphenol 
polymers. 


15 


FIGURE  5.  Infrared  Spectrum  of  Sample  Fraction  of  Compound  (I)- Biphenol  Polymer 


TABLE  2.  Results  of  Fractionation  of 
Compound  (I)-Biphenol  Polymer 


Weight  of 
Fraction,  g 

Percent 
of  Total 

Melting 
Range,  *C 

Molecular 

Weight 

2.2089 

56 

160-185 

18,100 

1.0195 

26 

160-170 

9.100 

0.2840 

7 

145-160 

5,800 

~0.1 

-3 

308-310 

1,600 

Polymerization  of  Hcxaphenyldichlorophospho- 
nitrile  Tetramer  with  Diamines 


Benzidine.  A  mixture  of  2.86  g  {0.004  mole)  of  compound  (I)  and 
0.73  g  (0.004  mole)  of  benzidine  was  heated  in  vacuum  to  300 #C  for  11 
hours.  During  the  heating  cycle,  83%  of  the  calculated  amount  of  hydro¬ 
gen  chloride  was  liberated.  The  product  was  a  brown  glass  which 
gradually  softened  over  the  range  205-250 *C  and  swelled  in  benzene  but 
did  not  dissolve. 

Anal.  Calc,  for  C.0H.AP.N. :  C.  69.9;  H.  4.8;  P.  15.1;  N,  10.2. 

Found:  C,  69.5;  H,  5.2;  P,  15.2;  N,  10.4;  Cl,  0. 

N.N1- Diphenylbenzidine.  A  mixture  of  3.36  g  (0.01  mole)  of  N,N'- 
diphenylbenzidine  and  7.13  g  (0.01  mole)  of  compound  (I)  was  heated  in 
an  evacuated  ampoule.  No  gas  evolution  occurred  until  about  275*C, 
and  the  reaction  became  vigorous  only  at  300  *C,  The  reaction  tempera¬ 
ture  was  held  at  360*C  for  6  hours.  About  90%  of  the  calculated  amount 
of  hydrogen  chloride  was  evolved. 

Anal.  Calc,  for  C, :  C,  73.8;  H.  4.9;  P.  12.7;  N,  8.6; 

'  60  48  4  6  ,  .  _ 

mol,  wt.  per  unit  976. 

Found:  C,  71.4;  H,  4.8;  P,  14.2;  N,  8.9; 

Cl,  0.5;  mol.  wt.  2090. 

Reaction  of  Compound  (1)  with  Trialkylphosphates 

A  mixture  of  2.14  g  (0.003  mole)  of  compound  (I)  and  0.60  g  (0.003 
mole)  of  triethylphosphate  was  heated  in  vacuum.  Gas  evolution  began 
around  150*C  and  the  mixture  changed  to  a  transparent  liquid.  Heating 
was  continued  for  10  hours;  the  temperature  was  gradually  raised  to 


300 *C.  The  volatile  product  (95%  of  theory)  was  identified  by  infrared 
analysis  as  ethyl  chloride.  The  residue,  a  light  yellow  resin,  melted 
at  125-145‘C. 


Anal.  Calc,  for  C, -H-.P.N.O 
-  38  35  5  4  4 

Found: 


C,  59.5;  H,  4.6;  P,  20.2;  N,  7.3; 
O.  8.4. 

C,  59.3;  H,  4.6;  P.  19.6;  N,  7.5; 
O.  8.8;  Cl,  0. 


Similar  results  were  obtained  with  both  tributylphosphate  and  tripropyl 
phosphate. 

Attempted  Reaction  of  Compound  (I) 
with  LiBr  in  Acetonitrile 


In  the  drybox,  a  1 -necked,  round-bottomed  flask,  equipped  with  a 
magnetic  stirrer  and  condenser,  was  loaded  with  10  g  (0.015  mole)  of 
compound  (1),  20  g  (O.n  mole)  of  LiBr  (dried  in  vacuum  at  200 *C  for 
18  hours),  and  500  ml  of  acetonitrile,  which  had  been  distilled  from 

p2o5. 

After  63  hours  at  the  reflux  temperature  of  CH3CN,  the  reaction 
mixture  was  filtered  free  of  solids.  Compound  (I),  9.8  g,  was  recovered 
unchanged  as  revealed  by  infrared  spectrum. 
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A  STUDY  OF  THE  THERMAL  DECOMPOSITION  OF  TRANS- 
2,6-DIAZIDOHEXAPHENY  LC  YCLOPHOSPHONITRILE 
TETRAMER  AND  FURTHER  STUDIES  ON  THE 
SOLVENT  POLYMERIZATION  OF  THE 
TETRAMER  WITH  BISPHOSPHINES 

by 

C.  M.  Sharts 
Chemistry  Division 


INTRODUCTION 

The  polymerization  of  trans-2,6-diazidocyclophosphonitrile  tetramer 
with  the  bisphosphines  1 ,4-bis(diphenylphosphino)butane  and  4,4'-bis- 
(diphenylphosphino)biphenyi  to  give  polymers  of  medium  chain  length 
was  reported  ea.lier.^ 

Also  reported  was  the  thermal  decomposition  of  trans-2,6-diazido- 
cyciophosphonitrile  tetramer.  This  thermal  decomposition  of  the  diazide 
tetramer  causes  considerable  concern  for  the  possibility  of  eventually 
obtaining  high  molecular  weight  polymers  from  the  trans-diazide  with 
bisphosphines.  A  decomposition  of  only  a  very  small  percentage  of  the 
starting  diazide  can  cause  premature  ending  of  polymer  chains  linked 
through  —  PsssN —  linkages. 

DISCUSSION 

An  investigation  has  been  made  of  the  melt-decomposition  of  trans- 
2.6-diazidocyclophosphonitrile  tetramer. ^  The  decomposition  was 
studied  with  a  thermal  balance  and  the  reaction  was  followed  by  means 
of  the  weight  loss  caused  by  nitrogen  evolution.  The  decomposition, 
shown  in  Table  3,  is  first  order  over  the  range  220-260*C.  The  activa¬ 
tion  energy  for  the  decomposition  was  37  Kcal/mole. 

Since  the  polymerizations  are  normally  carried  out  in  solution,  a  study 
was  made  of  the  decomposition  of  2.6-diazidocyclophosphonitrile  tetramer 
at  213*C  in  dibutylphthalate  solvent,  in  a  calibrated  vacuum  system. 

The  amount  of  data  collected  for  several  runs  was  voluminous.  Although 
final  calculations  to  the  three  significant  figures  justified  by  the  data 
have  not  been  completed,  preliminary  calculations  to  one  and  two  signi¬ 
ficant  figures  have  been  made  and  the  following  conclusions  reached: 


^NAVWEPS  Report  8226.  Quarterly  Report:  Foundational  Research 
Projects.  October-December  1964.  p.  16. 

^C.  M.  Sharts.  D.  R.  Gentry,  and  A.  Bilbo,  unpublished  results. 
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TABLE  3.  First  Order  Decomposition  of  trans-2,6 
Diazidocyclophosphonitrile  Tetramer 


Sample  size,  g 

Temperature,  *C 

iz  •  ”* 

K,  min 

0.314 

220 

1.6  X  10"2 

0.410 

240 

7.7  X  10‘2 

0.312 

240 

6.9  X  10'2 

0.144 

260 

28  x  10"2 

The  decomposition  is  first  order  in  dibutylphthalate  solution  with  a  half- 
life  at  213  #C  of  31  ±  2  min,  which  corresponds  to  a  rate  constant 

=  2.2  X  lO"^  min”*.  This  first-order  constant  at  213‘C  is  slightly 
greater  than  the  rate  constant  for  the  melt  decomposition  at  220 °C 

=  1.6  X  10"^  min”*).  The  slight  discrepancy  can  be  attributed 
to  the  difference  in  medium.  The  close  agreement  between  Kj^*^  for 
solution  decomposition  and  for  melt  decomposition  is  encouraging 

enough  to  permit  the  discussion  which  follows. 

It  seems  reasonable  to  assign  an  activation  energy  to  the  solution 
decomposition  of  diazide  identical  with  the  activation  energy  for  melt 
decomposition.  An  extrapolation  of  a  in  K  (K  =  rate  constant)  vs  l/T 
from  K|213  for  diazide  decomposition  with  an  assumed  activation  energy 
of  37  Kcal/mole  will  permit  an  estimation  of  the  rate  of  decomposition  of 
diazide  at  lower  temperatures.  The  results  of  such  a  plot  are  shown  in 
Table  4. 

The  last  column  of  Table  4  is  the  significant  one  for  consideration  of 
the  temperature  of  polymerization  of  trans-2,6-diazidocyclophospho- 
nitrile  tetramer.  It  represents  the  time  required  to  decompose  0.1% 
of  the  diazide  present.  For  the  chlorobenzene  solution  polymerization 
reported  earlier, ^  a  typical  time  for  96%  polymerization  was  on  the 
order  of  250  minutes.  At  the  132 *C  temperature,  the  time  for  0.1% 
decomposition  is  around  40  minutes.  In  other  words,  the  polymerization 
in  chlorobenzene  will  of  necessity  compete  with  an  amount  of  decompo¬ 
sition  which,  although  insignificant  for  normal  organic  reactions,  is 
very  important  for  polymerizations. 

Before  going  further,  it  is  important  to  recall  the  long  extrapolation 
that  has  been  made  to  get  the  times  for  0.1%  decomposition.  Table  4 


^NAVWEPS  Report  8226,  op  cit. 
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TABLE  4.  Extrapolated  Rate  Constants  and  Half-Lives  at 
Lower  Temperatures  for  Decomposition  of  2,6-Diazido- 
cyclophospbonitrilc  Tetramer  in  Dibutylphthalate 


Temperature, 

°C 

K  extrapolated, 

.  -1 
min 

Relative 

rate 

V* 

min 

S/IOOO 

min 

213 

2.2  X  10~2 

2.2  x  104 

31 

0.045 

202 

1.0  x  l(f2 

1.0  x  io4 

69 

0.10 

186 

2.2  X  10"3 

2.2  X  103 

310 

0.45 

175 

1.0  x  io-3 

1.0  X  JO3 

690 

1.0 

152 

2.2  X  10'4 

2.2  x  102 

3,100 

4.5 

137 

1.0  X  io‘4 

1.0  X  io2 

6,900 

10 

130 

2.2  X  l(f 5 

2.2  X  10 

31,000 

45 

118 

1.0  x  10"5 

1.0  X  10 

100 

111 

99 

2.2  X  10"6 

1.0  x  io-6 

2.2 

1.0 

310,000 

450 

1000 

can  be  wrong  by  an  order  of  magnitude.  Regardless  of  the  absolute 
accuracy  of  the  table,  the  results  clearly  show  that  lowering  the  tem¬ 
perature  will  be  extremely  advantageous  in  preventing  thermal  decom¬ 
position  of  the  diazidophosphonitrile  tetramer.  For  polymerizations 
whose  rate  of  reaction  is  such  that  3-4  hours  are  required,  it  would  be 
highly  desirable  for  temperatures  to  remain  below  115°C. 

A  number  of  polymerizations  were  carried  out  in  solvents  and  the 
results  of  the  runs  are  shown  in  Table  5.  Of  these,  only  Run  V  will  be 
given  detailed  treatment  in  the  Experimental  part  of  this  report. 

Runs  III  and  IV  are  noteworthy  in  that  polymers  of  molecular  weights 
9100  and  8900  were  obtained  for  the  diazide  tetramer  and  bisphosphine 
in  the  solvents  dibutylphthalate  and  xylene,  respectively.  Also  signifi¬ 
cant  in  Run  III  was  the  low  temperature  used,  110~115*C.  The  polymer 
was  white  when  isolated. 
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TABL.E  5.  Polymerization  of  Diazide  Tetramer  and  Bisphosphines  in  Solvents 
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The  polymer  was  white  until  heated  to  185"C.  About  0.9  g  was  ootained  — a  near -quantitative  yield. 
Xylene  worked  well.  Toluene  should  be  used,  or  possibly  heptane. 


The  most  promising  polymeric  material  of  Run  V  (sample  1076-42B- 
End)  was  subjected  to  thermal  gravimetric  analyses  in  both  helium  and 
air  (Figures  6  and  7)  and  found  to  be  stable  to  400 8C.  Decomposition 
was  slow  between  400*  and  450  °C  but  quite  rapid  between  450°  and  600  *C 
where  nearly  50%  of  the  weight  loss  occurred.  A  DTA  thermogram  (Fig¬ 
ure  8)  on  the  'ample  was  in  very  close  agreement  with  the  TGA  results. 

Thus  it  can  be  concluded  that  polymers  from  2,6-diazidocyclophospho- 
nitrile  tetramer  and  bisphosphines  remain  highly  promising  as  tempera¬ 
ture-stable  material.  Future  polymerization  studies  will  be  directed  to 
polymerizations  occurring  at  90- 100 *C.  Heptane,  toluene,  isopropyl  alco¬ 
hol,  and  a  suitable  ester  will  be  investigated  as  solvents  for  polymerization. 

EXPERIMENTAL 

Low  Temperature  Polymerization  of  trans-2t6-Diazido- 
cyclophosphonitrile  Tetramer  and  4,4t-bis(diphenylphos- 
phino)biphenyl  in  Dibutylphthalate  Solvent,  Run  V 

Diazide  tetramer  (1.0885  g,  1.498  mmole),  bisphosphine  (0.7827  g, 

1.498  mmole),  and  60  ml  of  dibutylphthalate  were  placed  in  a  flask  which 
was  then  evacuated  to  degas  solvent  and  reactants.  After  extended  (1  hour) 
degassing,  the  contents  of  the  flask  was  heated  slowly.  Nitrogen  evolution 
was  observed  at  95 *C.  The  temperature  was  maintained  overnight  between 
95°  and  110*C  and  during  the  night,  white  polymer  precipitated.  The  sol¬ 
vent  was  slightly  discolored  to  a  pale  yellow.  About  0.11  g  of  diazide  was 
added  to  the  flask.  Dibutylphthalate  was  then  distilled  off  under  reduced 
pressure.  Care  was  taken  not  to  exceed  160*C  pot  temperature  at  any 
time.  The  polymer  residue  should  have  been  all  azide-terminated  at  this 
time.  The  reaction  flask  was  maintained  at  160 *C  for  several  hours  to 
assure  complete  reaction  of  all  phosphine  end  groups.  Most  of  the  poly¬ 
meric  solid  dissolved  when  it  was  treated  with  50  ml  of  hot  benzene.  The 
benzene-insoluble  polymer,  amounting  to  about  0.2  g,  was  labelled  "Sample 
A"  and  set  aside  to  dry.  Upon  the  addition  of  100  ml  of  heptane  to  the 
benzene  solution,  polymer  precipitated  which  was  labelled  "Sample  B." 

The  material  in  the  benzene-heptane  solution  was  recovered  as  an  oil  and 
labelled  "Sample  C."  Samples  A,  B,  and  C  were  each  dissolved  in  o-di- 
chlorobenzene  and  0.1  g  of  triphenylphosphine  was  added  to  each.  The 
endcapping  reaction  was  carried  out  overnight  and  the  heptane  was  added 
to  the  hot  o-dichlorobenzene  solutions.  Endcapped  polymers  were  obtained 
which  had  the  final  designations  1076-42A  End,  1076-42B  End,  and  1076-42C 
End.  Analyses  and  molecular  weight  determinations  of  these  samples  are 
pending. 


The  thermal  gravimetric  analyses  (TGA)  were  made  on  a  Stanton 
Model  H.T.  thermal  balance.  The  differential  thermal  analyses  (DTA) 
were  obtained  on  a  Du  Pont  900  DTA. 
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FIGURE  8.  Differential  Thermal  Analysis  of  Sample  1076-42B-End  in  Nitrog 


INFRARED  ATOMIC  SPECTRA 


THE  FIRST  SPECTRUM  OF  MERCURY  IN  THE 
3.2-  TO  4.0-MICRON  REGION 

by 

C.  J.  Humphreys  and  E.  Paul,  Jr. 
Research  Department 


INTRODUCTION 

The  spectrum  of  neutral  mercury,  Hg  I,  is  one  of  the  best-described 
and  analyzed  among  atomic  emission  spectra.  It  was  also  one  of  the 
first  to  be  studied,  being  included  in  the  investigations  of  Kayser  and 
Runge  and  of  Rydberg  before  the  turn  of  the  century.  The  results  of 
these  early  investigations  are  included  and  referenced  in  the  monograph 
by  Paschen  and  Gotze,*  as  are  also  Paschen's  own  observations  in  the 
extra- photographic  infrared  region. 

One  of  the  authors  of  this  section  published  a  paper  in  1953, ^  reporting 
and  discussing  observations  of  Hg  1  made  somewhat  earlier  at  the  National 
Bureau  of  Standards  and  covering  the  region  between  1.3  and  2.0  microns 
accessible  to  observation  with  an  ambient-temperature  lead  sulfide  detec¬ 
tor.  References  to  the  then  available  literature  on  the  analysis  of  Hg  I 
were  included.  Most  of  the  lines  reported  in  that  paper  had  already  been 
observed  during  the  Paschen  era.  The  principal  new  experimental  con¬ 
tribution  was  the  resolution  of  the  satellite  lines  of  the  6d*D—  5f^F* 
multiplet. 

The  mercury  level  system  is  similar  to  that  of  the  other  homolo¬ 
gous  elements,  zinc  and  cadmium,  but.  as  might  be  expected  in 
the  instance  of  an  element  of  so  great  an  atomic  mass,  there  is 
considerable  departure  from  L S  coupling.  Figure  1  is  a  Grotrian 
diagram  showing  the  portions  of  the  energy  level  scheme  that  are  required 
for  discussion  of  the  material  obtained  in  the  experiments  herein  reported. 
Further  discussion  will  be  given  of  the  observed  transitions  illustrated 
in  the  figure.  The  normal  configuration  d*®s^  accounts  for  the  lowest 
term  *Sq.  Excitation  of  an  electron  from  s^  to  a  succession  of  s.  p.  d. 


*F.  Paschen  and  R,  Gdtse,  Seriengesetze  der  Linienspektren,  Berlin: 
Springer.  1922. 

7 

C.  J.  Humphreys,  J.  Opt.  Soc.  Am..  Vol,  43,  p.  1027  (1953). 
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FIGURE  1.  Partial  Energy* I  vel  Diagram  of  Hg  I.  Circle*  represent 
odd  levels  and  horisontal  lines  even  levels.  Some  of  the  more  intense 
transitions  beyond  3  microns  are  indicated. 


f,  or  g  states  gives  rise  to  long  series  of  triplet  and  singlet  levels,  all 
of  which  converge  to  the  d^s^SQl  normal  state  of  the  mercury  ion. 

This  simple  pattern  is  broken  by  the  appearance  of  levels  from  the  con¬ 
figuration  resulting  from  binding  an  electron  to  the  d^s^  configuration 
of  the  ion.  The  d^s^  configuration  accounts  for  ar.  inverted  term  in 
the  ion  spectrum.  Owing  to  extremely  large  spin-orbit  interaction  in 
the  ion,  the  ion  doublet  separation  is  very  large,  amounting  to  15038 
cm"*.  All  of  the  Hg  I  levels  in  this  family  that  lead  to  high-energy  tran¬ 
sitions  arise  from  the  configuration  d^s^  •  p  and  converge  to  the  lower 
doublet  limit  ^D2±.  Such  known  levels,  conventionally  designated  by  a 

prime,  are  6p,3P2»  6p'3Dj,  6p,3F4«  6p'^D2,  and  6p^PJ.  Levels  asso¬ 
ciated  with  the  upper  level  of  the  inverted  ^D,  although  known,  are  not 
of  significance  in  connection  with  the  material  introduced  into  this  report. 
The  levels  of  the  d^s^(^D2|)np  family  are  intermingled  with  those  built 
upon  d^s(^S).  As  a  result,  the  mercury  spectrum  is  characterized  by 
numerous  perturbations  between  levels  of  the  same  parity  and  j-value. 

The  most  conspicuous  example  is  the  perturbation  of  Sf^F^  by  bp'^F^, 
which  leads  to  extremely  irregular  intervals  in  the  5f3F  multiplet.  This 
situation  was  discussed  in  the  earlier  paper,3  but  it  is  also  significant 
for  the  identification  and  evaluation  of  the  complex  level  5gG,  which  is 
one  of  the  items  with  which  this  report  is  concerned. 

Recent  observations  of  the  mercury  spectrum  have  been  to  a  great 
extent  features  of  various  laser  programs.  The  observations  described 
and  discussed  here  were  stimulated  by  receipt  of  an  advance  copy  of  a 
paper  by  Bockasten,  Garavaglia,  Lengyel.  and  Lundholm  that  appeared 
subsequently  in  the  Journal  of  the  Optical  Society  of  America.*  That 
paper  reported  laser  lines  in  Hg  I  in  the  interval  between  1  and  2.5 
microns.  Discussion  was  included,  however,  of  several  laser  transitions 
reported  by  other  observers  at  wavelengths  greater  than  2.5  microns. 
Included  among  these  was  a  feature  at  3.932  p  first  observed  in  spontane¬ 
ous  emission  by  Paschen^  and  classified  by  him  ambiguously  as 
bp,3P^-bd3D3  or  5fF*  — 5gG.  On  the  basis  of  this  background  material, 
it  appeared  desirable  to  try  to  determine  the  origin  of  the  feature  by 
observing  it  with  higher  resolution  than  had  been  employed  heretofore. 
Preliminary  results  of  this  study  were  reported  at  the  Philadelphia 
meeting  of  the  Optical  Society  of  America,  October  1965.^  Subsequent 
to  submission  of  the  abstract,  the  problem  was  discussed  with  Dr. 


3Ibid. 

K.  Bockasten.  M.  Garavaglia,  B.  A.  Lengyel,  and  T.  Lundholm, 

J.  Opt.  Soc.  Am.,  Vol.  55,  p.  1051  (1965). 

^F.  Paschen,  Ann.  Physlk,  Vol.  6,  p.  47  (1930). 

^C.  J.  Humphreys  and  E.  Paul.  Jr.,  J.  Opt.  Soc.  Am.,  Vol.  55, 
p.  1572  A  (1965). 
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Bockasten  in  Sweden  and  with  k.  C.  M.  Learner  of  the  Imperial  College, 
South  Kensington.  Dr.  Learner  was  interested  in  the  possibility  f 
observing  the  transition  bp’^F^—  5gG,  assuming,  of  course,  that  the 
5gG  complex  had  been  found.  The  expected  location  was  at  about  3.522  p. 
In  brief,  the  investigation  demonstrated  that  the  ambiguously  classified 
feature  at  3.932  p  consisted  of  four  resolved  components,  three  of  which 
represented  a  portion  of  the  5fF  — 5gG  multiplet,  the  fourth  and  strongest 
being  bp'^P^  —  bd^Dj.  The  transition  bp'^F^— 5gG  was  also  observed  as 
a  moderately  intense  line.  Significant  details  are  reported  below. 

EXPERIMENTS  AND  ANALYSIS 

Observations  of  the  infrared  spectrum  of  mercury  were  made  by 
means  of  a  *^®Hg  microwave-excited  electrodeless  source  and  the  same 
grating  spectrometer  (mounting  a  75G0-!ines-per-inch  grating)  ■  -d  in 
all  recently  reported  experiments.  Although  the  problems  of  ediate 
interest  involved  only  the  region  between  3.5  and  4.0  p,  the  entii  r-mge 
from  1.8  to  a  little  beyond  4.0  p  was  covered.  A  section  of  chart  with 
identified  features  covering  the  region  between  3.9  and  4.0  p  is  repro¬ 
duced  as  Figure  2.  The  group  of  four  lines  at  the  short  wavelength  end 
of  this  figure  is  identified  with  the  ambiguously  classified  feature  at 
3.93  p.  The  wavelength  of  the  line  6p,3p2  —  6d3D3  originating  in  I 

is  calculable  from  levels  derived  from  interferometric  measurements.^’® 
Measurement  of  the  positions  of  the  other  three  lines  of  this  group, 
together  with  that  of  the  transition  identified  as  5f3F£— -5gG,  permits 
evaluation  of  the  complex  of  energy  levels  described  as  5gG.  Agreement 
between  determinations  of  this  value,  79783.85  cm**,  from  the  four 
available  transitions  is  close  enough  to  estimate  the  precision  at  ±0.02 
cm**.  It  is  noted  at  this  point  that  the  separation,  amounting  to  about 
50  cm’1,  of  the  level  5f3F}  from  the  group  of  closely  spaced  levels, 
5f3F|.  5f3FJ,  and  5f*Fj,  is  anomolcus,  probably  being  caused  by  a 
perturbation  between  5f  3F4  and  6p,3FJ.  Inspection  of  higher- series 
f-levels  reveals  no  s«ch  anomaly.  Keeping  in  mind  the  fact  that  the 
f-level*  should  show  very  small  intervals,  it  is  not  unexpected  that  the 
5gG  levels  should  be  completely  unresolved  within  tht  limits  of  available 
instrumental  capabilities. 

The  levels  of  *^®Hg  I  that  are  involved  either  directly  or  indirectly 
in  the  determination  of  5gG  are  compiled  in  Table  1.  A  description  of 
the  spectrum  between  3.2  p  and  the  upper  wavelength  limit  of  the  obser¬ 
vations  is  compiled  and  presented  as  Table  2.  Unfortunately,  not  all 

^K.  Burns  and  K.  B.  Adams,  J.  Opt.  Soc,  Am..  Voi.  42,  p,  56  (1952). 

8  ” 

C.  J.  Humphreys  and  E.  v’aul,  Jr.,  J.  Physique,  Vol.  19.  p.  424 

(1958);  Trans.  I.  A.  U.,  Vol.  1 1  A,  p.  07  (1962). 
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TABLE  1.  Selected  Levels  in  Hg  I 

Designation 

-1* 

Level  value,  cm 

v3p; 

68886.445 

6d*D2 

71333.334 

6d3DJ 

71336.299 

6d3D2 

71396.368 

6d3E3 

71431.466 

5f3F- 

77237.016 

3 

5f  F' 

77239.344 

5f1F* 

77241.656 

5f3F  * 

4 

77287.027 

5gG 

79783.85  ±  0.02 

The  value  of  5gG  is  based  on  current  observa¬ 
tions;  all  others  are  derived  from  interferometric 
measurements  by  Burns  and  Adams  and  by 
Humpnreys  and  Paul. 
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TABLE  2.  Description  of  Hg  land  Hg  I,  Natural 
(3.2-4  Micron  Region) 


Estimated 

r~ - - - - 1 

X, 

• 

c r, 

vac.,  cm"* 

Transition 

intensity 

air,  A 

600C 

32135.674 

3110.958 

7p*P1‘-8»1S0 

250 

34724,16 

2879.05 

8p3P‘-8d3D3 

30 

34872.98 

2866.77 

8p3P*-8d3D2 

5000 

34927.78 

2862.27 

88  3SJ  —8p3P® 

200 

*  35217.5 

2838.7 

6p,3F®  —  5gG 

150 

t  35736.18 

2797.522 

8p1P1°-8d1D2 

7000 

36303.07 

2753.84 

7p3P'-8s3S1 

't  37286.6 

2681.2 

7diD2  -6*^* 

200 

< 

t  37290.7 

W 

2680.9 

160 

t  37308.8 

2679.6 

7d‘D2-6f3F- 

100 

37475.98 

2667.65 

7p1P®-8s3S1 

300 

t 37596.5 

2659.1 

7d3D1-6f3F® 

"t  37905.8 

2637.4 

7d3D2-6f1F® 

80 

X 

t  37910.1 

v. 

2637.1 

7d3D2-6f3F® 

200 

t  37928.8 

2635.8 

7d3D2-6f3F* 

800 

t  37175.0 

2618.8 

7d3D  -6f3F® 

3  4 

600 

*39254.67 

2546.77 

5f3F®  — 5gG 

5000 

39281.849 

2545.011 

6p'3P®-6d3D3 

1500 

*39290.11 

2544.48 

5f3F“-5gG 

(Footnote  at  end  of  table) 


TABLE  2  (Contd.) 


Estimated 

intensity 

X. 

© 

air,  A 

o-, 

vac.,  cm" * 

Transition 

1000 

*  39324.70 

2542.24 

5f*  F®  —  5gG 

300 

39831.157 

2509.913 

6p'3P£-6d3D2 

900 

, 

t  39897.029 

2505.769 

8s3S1-8p3P1° 

700 

i 

*40039.50 

2496.85 

5f3F“-5gG 

200 

t  4021 5.1 85 

2485.945 

8sV8P3po 

400 

t 40661.242 

2458.674 

8s1S0-8p*PJ 

50 

40857.246 

2446.879 

6p*3P*~6d1D2 

Note:  Unmarked  wavelengths  are  calculated  from  levels. 

*  Measured  wavelengths.  source, 

t  Calculated  from  levels  of  natural  Hg. 


observed  transitions  can  be  evaluated  from  available  energy  levels 
obtained  from  observations  on  ^®Hg  I.  In  such  instances  calculated 
wavelengths  are  those  of  natural  mercury  obtained  from  levels  listed  in 
"Atomic  Energy  Levels. Apology  is  made  for  publishing  this  mixed 
list,  but  it  was  the  best  that  could  be  done  since  ^®Hg  was  used  as  the 
source. 

CONCLUSION 

The  interpretation  of  a  feature  (namely,  \3.93  p)  obtained  in  spontane¬ 
ous  emission  and  also  observed  as  a  laser  transition  leaves  an  obvious 
question  open.  If  laser  emission  were  observed  with  sufficient  resolving 
power,  would  both  categories  of  transitions  observed  in  spontaneous 
emission  also  appear,  or  to  what  extent  would  the  entire  spontaneously 
emitted  spectrum  be  reproduced?  To  answer  these  questions  it  is  highly 
desirable  that  further  observations  with  lasers  coupled  with  high- reso¬ 
lution  equipment  be  undertaken. 


q 

C.  E.  Moore,  "Atomic  Energy  Levels,"  Circular  of  the  National 
Bureau  of  Standards  467,  Vol,  III,  p.  192  (1958). 
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NOTE  ON  CESIUM  I 

by 

C.  J.  Humphreys  and  E.  Paul,  Jr. 
Research  Department 


The  first  spectrum  of  cesium  is  noteworthy  because  of  the  relatively 
small  interval  separating  the  and  terms  comprising  the  first 
excited  levels  above  the  normal  '^S  state.  Owing  to  the  relative  sim¬ 
plicity  of  the  spectrum  with  the  energy  distributed  over  a  small  number 
of  transitions,  combinations  of  such  low  terms  are  expected  to  be  intense. 
Accordingly,  the  multiple!  of  three  lines  comprising  the  permitted  tran¬ 
sitions  between  the  levels  of  the  listed  terms,  classically  described  as 
the  diffuse  doublet,  is  not  only  located  relatively  far  out  in  the  infrared 
but  is  probably  the  most  intense  group  of  atomic  emission  lines  to  be 
found  in  the  given  region.  Based  on  the  best  available  level  values, 
which  require  further  discussion,  the  wavelengths  of  the  three  lines 
are  as  follows: 


6p2p4 
&p  2pU 


e  ,2_ 

-5d  D  i 

1  2 

30103.355  A, 

air 

2 

0 

—  5d  D  i 

34900.137  A, 

air 

-5d2D  i 

A  2 

36131.089  A, 

air 

The  cesium  levels  are  characterized  by  a  hyperfine  structure  splitting, 
but  this  splitting  is  too  small  to  be  observed  by  direct  observation  except 
in  the  instances  of  ^SqI  and  ^Pq|.. 

Because  of  its  intensity,  location,  and  essentially  single  character, 
the  strongest  line  of  the  group,  X.3.49  n,  has  appeared  to  be  a  good  candi¬ 
date  for  precise  evaluation  as  a  wavelength  standard.  During  the  summer 
of  1964,  when  the  infrared  spectrum  of  f36Xe  was  being  measured  inter- 
ferometrically  in  preparation  for  the  Hamburg  meeting  of  the  International 
Astronomical  Union,  an  attempt  was  made  to  intercompare  Cs,  X.3.49  n* 
with  several  *  3^Xe  lines  in  the  same  region.  At  that  time  the  prelimi¬ 
nary  wavelength  of  Cs,  X3.49  p,  was  calculated  from  levels  published 
by  H.  Kleinian.*®  Several  recordings  were  made  using  an  etalon  of  30 
mm.  Because  it  was  not  found  possible  to  establish  the  whole  order  of 
interference  for  the  Cs  patterns,  the  attempt  failed,  suggesting  that  the 
provisional  wavelength  was  in  error.  Subsequently,  Eriksson,  Johansson, 


H,  Kleiman,  J.  Opt.  Soc.  Am., 


Vol.  52,  p.  441  (1962). 
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and  Norlen*  *  remeasured  the  cesium  spectrum  and  obtained  revised 
values  of  the  6p^P  levels,  indicating  that  a  portion  of  Kleiman's  level 
determinations  had  been  in  error.  On  the  basis  of  this  new  information, 
the  wavelength  of  X.3,49  p  was  recalculated.  The  old  records  were  then 
re-examined.  It  was  found  that  intercomparison  of  three  *  ^Xe  iine3i 
namely,  35070.2520,  33666.6991,  and  32739.2788  A,  led  to^agreement 
with  the  calculated  value  for  the  Cs  wavelength  34900.137  A,  within  one 
in  the  thousandths  place.  Actually,  the  reductions  led  to  the  following 
uncorrected  values  listed  opposite  the  appropriate  ^^Xe  line: 


Cs 


35070.2520  A 
33666.6991  A 
32739.2788  A 


34900.1380  A 
34900.1366  A 
34900.1365  A 


This  agreement  constitutes  a  verification  of  the  correctness  of  the 
observations  by  Eriksson,  Johansson,  and  Norlen.  To  complete  the 
record,  the  following  tabulation  is  given  of  the  levels  used  in  the  calcu¬ 
lation  of  diffuse  doublet  lines  in  Cs  I: 


&P2p0l 

11178.2697 

Center  of  gravity,  E,  J,  and  N 

6p2p;i 

11732.3079 

Center  of  gravity,  E,  J,  and  N 

5d2D  , 

1  2 

14499.253 

Center  of  gravity,  Kleiman 

5dZD  , 

^2 

14596.845 

Center  of  gravity,  Kleiman 

Additional  intercomparisons  making  use  of  longer  interference  paths 
are  desirable.  These  should  include  the  other  cesium  lines  of  the  p-d 
multiplet  as  well  as  still  others  in  the  same  region. 


11  <* 

K.  B.  Eriksson,  I.  Johansson,  and  G.  Norlen,  Arkiv  for  Fysik, 

Vol,  28,  p.  233  (1964).  “  '  — 
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PLASMA  PHYSICS 


PRODUCTION  OF  HIGH  POWER  PULSES 
BY  MAGNETIC  FIELD  COMPRESSION 

by 

R.  L.  Conger,  J.  H.  Johnson,  L.  T.  Long,  and  J.  A.  Parks 
Electricity  and  Magnetism  Division 


INTRODUCTION 

The  object  of  this  research  is  to  investigate  the  possibility  of  con¬ 
structing  a  small  device  that,  with  energy  provided  from  a  chemical 
explosive,  could  produce  10^  joule  pulses  of  electricity  at  10’  watts. 

A  transient  field  of  the  order  of  100  kG,  lasting  for  a  few  hundred 
microseconds,  can  be  produced  by  discharging  a  capacitor  bank  into  a 
strongly  constructed  solenoid  of  a  few  turns.  Higher  fields  cannot  be 
obtained  because  the  forces  produced  by  the  magnetic  field  destroy  the 
solenoid.  However,  fields  of  over  a  million  gauss  can  be  produced  by 
using  high  explosives1  to  implode  a  one-turn  coil  in  the  form  of  a  cylin¬ 
der.  If  the  implosion  takes  place  in  a  time  small  compared  to  the  decay 
time  of  current  in  the  cylinder,  the  total  flux  will  be  conserved  and  the 
magnetic  field  will  be  inversely  proportional  to  the  decreasing  volume 
of  the  cylinder.  Cold  rolled  copper  and  stainless  steel  have  a  high 
enough  conductivity  for  field  compressors  of  this  type.  Although  the 
cylinders  are  usually  circular,  they  can  be  of  other  shapes. 

Since  the  energy  in  the  magnetic  field  is  given  by  the  integral^ 

l  =  ijt iH2dV 

and  the  field  is  roughly  inversely  proportional  to  the  volume,  the  energy 
in  the  field  compressor  is  also  inversely  proportional  to  the  volume  and 
therefore  is  increased  in  the  implosion.  This  increase  comes  from  the 
chemical  energy  of  the  explosive.  If  the  initial  field  before  the  implosion 


*R.  S.  Caird,  W.  B.  Gam,  D.  B.  Thomson,  and  C.  M.  Fowler,  J.  Appl. 
Fhys.,  Vol.  35,  p.  781  (1964). 

I 

J.  A.  Stratton,  Electromagnetic  Theory,  International  Series  in 
Pure  and  Applied  Physics.  New  York:  McGraw-Hill,  1941,  p.  124. 
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is  sufficiently  high,  10  percent  or  more  of  the  energy  in  the  explosive 
can  be  transferred  to  the  magnetic  field  energy.  At  the  limit  of  the 
implosion,  the  magnetic  pressure  exceeds  the  pressure  produced  by  the 
explosive  and  the  field  starts  to  expand  again.  Since  the  current  in  the 
compressor  is  proportional  to  the  magnetic  field,  this  current  is  also 
inversely  proportional  to  the  volume  of  the  compressor. 

Although  the  magnetic  field  compressor  converts  the  chemical  energy 
of  the  explosive  to  electrical  energy  with  high  efficiency,  getting  this 
energy  out  can  be  difficult.  The  conductance  of  the  field  compressor 
plus  any  constant  load  applied  during  the  compression  must  be  sufficiently 
high  so  that  the  decay  time  of  the  magnetic  field  is  much  longer  than  the 
time  of  the  implosion.  More  efficient  conversion  of  the  chemical  energy 
is  obtained  if  the  field  is  compressed  with  little  or  no  energy  removed, 
and  then,  after  the  field  is  compressed,  the  current  is  switched  to  the 
load.  There  are  several  possible  ways  of  accomplishing  this.  One 
method  is  to  place  a  strip  of  tungsten  in  the  circuit  of  the  field  compres¬ 
sor,  or  some  similar  material  which  has  a  high  conductivity  when  cold 
but  has  a  small  enough  mass  so  that  it  is  heated  to  a  high  temperature 
by  the  high  current  in  the  compressor  and  strip.  As  it  heats  up,  the 
resistance  of  the  strip  increases,  thus  switching  the  current  to  the  load. 

Since  the  total  magnetic  flux  in  this  device  is  conserved,  it  is  not 
possible  to  cascade  field  compressors  to  produce  larger  and  larger  fields; 
however,  it  is  possible  to  obtain  energy  from  each  stage  of  a  cascaded 
field  compressor  as  the  field  is  allowed  to  expand  into  the  next  stage. 

The  first  step  in  developing  such  a  cascaded  compressor  would  be  to 
develop  a  single  stage  in  which  energy  could  be  removed  efficiently. 


THEORY  OF  MAGNETIC  FIELD  COMPRESSION 


Magnetic  field  compression  calculations  can  be  based  on  the  integral 
form  of  Maxwell's  equation; 
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where  E  is  the  electric  field,  B  is  the  magnetic  field,  t  is  the  time,  and 
cgs  units  are  used. 


From  a  vector  identity  and  Ohm's  law, 
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where  R  is  the  resistance  of  the  cylinder  in  which  the  field  is  being 
compressed  and  i  is  the  current. 

The  field  compressor  could  be  constructed  as  shown  in  Figure  l. 
The  explosive  would  cause  the  short  to  advance  to  the  right  at  the  deto¬ 
nation  velocity  v. 


SHORT  PLATES  PLATE  RESISTANCE 


FIGURE  1.  Flat  Magnetic  Field  Compressor 


For  this  geometry  equations  (1)  and  (2)  give  the  differential  equation 

+iBf  =  iR*‘°8 

If  edge  corrections  are  neglected 


and  thus 


SwRdtXlO®  dX  dB 
2ira(XQ  -  X)  *  XQ  -  X  ‘  B 


The  advancing  shock  wave  from  the  explosive  drives  the  conductors 
of  Figure  1  together  at  a  velocity  v,  so  that  X  =  vt.  Substitution  of  this 
value  for  X  in  equation  (5),  followed  by  integration,  gives 


B  = 


where 


P  = 


5wR  X  IQ8 
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(7) 


If  R  =  0,  then  B  is  inversely  proportional  to  1  -  X/Xq,  and  the  field 
is  inversely  proportional  to  the  volume  between  the  plates;  but  for  R 
finite  some  of  the  magnetic  field  leaks  through  the  conductor.  If  0  >  1, 
no  field  compression  takes  place. 


Since  B  is  proportional  to  i,  equation  (6)  car.  be  written  in  terms  of 
current  to  give 


i 


(8) 


For  p  »  1  and  X  «  Xq,  the  logarithm  of  equation  (8)  gives 


(9) 


or 


■|Jvt/X. 


a  =  iQe 


(10) 


The  inductance  of  the  field  compressor  shown  in  Figure  1  is  given  by 


L  = 


2wa(XQ  -  X)  x  10 
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if  edge  connections  are  neglected.  From  equations  (7)  and  (11) 


RX_ 


(12) 
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where  Lq  is  the  value  of  L  when  X  =  0.  For  X  «  Xq  and  by  the  use  of 
equation  (12),  equation  (10)  can  be  written  as 


i  2 


-Rt/L 

e 


(13) 


which  is  the  usual  equation  for  the  decay  of  current  in  an  inductance. 
Field  compression  must  take  place  in  a  time  that  is  small  compared  to 
the  L/R  time  constant  of  the  field  compressor. 

When  p  is  expressed  in  terms  of  equation  (12),  equation  (6)  can  be 
applied  to  geometries  other  than  the  one  shown  in  Figure  1.  In  particu¬ 
lar,  it  can  be  applied  to  a  helical  coil  with  a  metal  tube  mounted  axially 
inside  the  coil.  In  a  field  compressor  of  this  form,  the  field  is  com¬ 
pressed  between  the  coil  and  the  tube  as  the  tube  is  expanded  explosively 
starting  at  one  end.  Since  this  type  of  field  compressor  can  have  a  much 
higher  inductance  than  the  one  shown  in  Figure  1,  the  impedance  can 
also  be  much  higher. 

When  the  resistance  R  of  the  field  compressor  is  constant,  the  energy 
transferred  to  it  is  given  by  the  integral 


i 


This  integral  can  be  maximized  by  differentiating  with  respect  to  R  and 
setting  the  result  equal  to  0.  If  this  is  done  it  is  found  that  W  is  a  maxi 
mum  when  p  ~  For  this  value  of  P,  equation  (14)  becomes 


0  0 


in  n 


where  n  is  the  compression  ratio  or  Xq/(Xq  -  X). 
From  equation  (12).  with  £  =  | 


V 
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Substitution  of  (16)  into  equation  (IS)  gives 

C  -  ii02LQlnn 


(15) 


(16) 
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Since  the  initial  energy  before  compression  is  }iQ^LQ,  the  energy  increase 
ratio  is  inn.  However  the  term  inn  does  not  increase  very  rapidly  with 
n.  More  efficient  transfer  of  energy  from  the  explosive  io  the  load  ca  i 
be  obtained  by  making  the  load  resistor  of  a  sheet  of  tungsten  or  tantalum 
thin  enough  to  heat  up  during  the  compression,  thereby  increasing  the 
resistance  and  transferring  energy  to  the  load. 

When  R  is  not  constant,  equation  {14)  cannot  be  evaluated  analytically. 
Two  methods  were  used  to  evaluate  this  equation  with  R  as  a  variable. 

One  method  was  the  use  of  an  analog  computer;  the  other  was  a  desk 
calculation  that  consisted  of  dividing  the  compression  into  a  number  of 
small  increments  and  letting  R  be  a  constant  in  each  increment.  The 
successive  values  of  R  were  calculated  from  the  energy  transferred  to 
R  in  the  preceding  increment.  Both  methods  gave  the  same  result,  with 
the  desk  calculation  serving  as  a  check  or.  *he  co  r  j-uter  calculation. 

Only  the  computer  calculations  will  be  discussed  in  this  report,  since 
they  are  probably  the  more  accurate  of  the  two. 

It  is  desired  to  maximize  the  energy  delivered  to  the  load.  The  equa¬ 
tion  for  this  energy  is 


l  * 


{!€) 


where  i  and  R  are  variables.  The  term  R  has  an  upper  limit  determined 
by  the  melting  point  of  the  material  used  for  the  load.  For  tungsten 
Rmax  =  12Rq,  and  for  tantalum  Rmax  3  9Ro.  where  R(j  is  the  initial 
resistance  of  the  load.  The  compression  ratio  determines  tj.  which, 
for  the  computer  calculations,  is  0.9$  Xq/v. 


The  variables  i  and  R  are  determined  by  the  pair  of  differential 
equations 


and 


? 

i  R 


C  .  dR 
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(19) 


(20) 


Equation  (19)  is  derived  from  equation  (3),  and  equation  (20)  represents 
the  change  in  R  caused  by  its  heating.  Here  Xq.  v,  L.  a.  and  Rq  are 
constants,  C  is  the  heat  capacity  of  the  load,  and  a  is  the  temperature 
coefficient  of  resistance  m  the  equation 
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R  =  RQU  +  qT) 


where  T  is  temperature. 

For  these  calculations,  the  constants  had  the  following  values: 

XQ  =  25  cm 

v  =  7  x  105  cm/sec 
-3 

L  *  4.4  x  1 0  h/  sec 
a  =  0.0045  ohm/*C  for  tungsten 
a  =  0.0035  ohm/*C  for  tantalum 
t,  =  3.392  x  10-5  sec 


=  6x10  amp 


The  computer  was  used  to  und  a  value  of  C/aRp  for  a  given  Rq  such 
that  R  =  Rmax  at  t  =  -.392  x  !0"5  sec.  When  this  constraint  was  met. 


the  value  of  E  for  that  value  of  C/aR(j  was  recorded.  Figure  2  presents 
£  as  a  function  of  Rq/L  for  the  case  of  Rmax  ~  12Rq,  and  E  as  a  function 
of  Rq/L  for  *  9Rq  i*  presented  in  Figure  3. 


These  figures  show  that  E  is  a  maximum  when 


ar  0.1 


From  equations  (18)  and  (20) 


CAR 


For  tungsten,  with  3^**  s  12Rq.  AR  *  11  Rq,  and 


For  tantalum,  with  R^^  ~  9  Rq,  AR  =  8Rq,  and 
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The  energy  £  in  equations  (24)  and  (25)  is  the  maximurr  value  of  the 
curves  of  Figures  2  and  3,  respectively.  From  the  values  of  Rq  from 
equation  (22),  and  C  from  equation  (24)  or  (25),  the  required  shape  of 
the  load  can  be  determined. 

Figure  4  shows  the  computer  mechanization  diagram  for  this  problem. 

MAGNETIC  FIELD  COMPRESSION  TESTS 

Figure  5  shows  the  magnetic  field  compressor  before  the  0.25  in. 
thick  penta-erythritol  tetranitrate  (PETN)  explosive  was  glued  to  the  top 
of  the  upper  copper  plate,  which  is  0.75  in.  above  the  lower  plate. 

A  steel  plate  1  in.  thick  above  the  explosive  partially  confines  the  explo¬ 
sive  and  therefore  increases  the  detonation  pressure  on  the  upper  copper 
plate.  The  three  RG-9B/U  coaxial  cables  attached  at  the  right  are  con¬ 
nected  to  a  480-|if  capacitor  bank  30  ft  away,  which  produces  the  initial 
field.  An  exploding-bridge-wire  (EBW)  type  detonator  at  the  top  of  the 
wooden  assembly  on  the  right  initiates  a  17.2  cm  length  of  explosive 
cord,  which  in  turn  initiates  the  PETN  explosive  after  a  36.5-psec  delay. 
This  delay  allows  the  magnetic  field  in  the  compressor  to  reach  its  maxi¬ 
mum  initial  value  before  the  detonation  wave  progressing  down  the  ex¬ 
plosive  forces  the  upper  copper  plate  against  the  lower  one.  The  wedge 
on  the  lower  plate  at  the  right  end  of  the  compressor  allows  contact  to 
be  made  between  the  upper  and  lower  copper  plates  when  the  detonation 
wave  has  traveled  5.5  cm.  Both  plates  are  15  in.  long,  1.5  in.  wide, 
and  0.070  in.  thick.  As  the  compression  progresses,  current  is  forced 
through  the  load,  which  is  a  0.007  in.  tantalum  strip  at  the  left  side 
behind  the  resistors.  The  tantalum  strip  is  placed  at  an  angle  of  30  deg 
so  that  there  will  be  less  dead  space  at  the  end  of  the  compression  when 
the  detonation  wave  reaches  the  left  end  of  the  upper  plate.  The  steel 
block  on  which  the  lower  plate  is  mounted  is  sufficiently  large  to  with¬ 
stand  the  force  of  the  explosion.  This  block  becomes  dented  and  the 
lower  copper  strip  is  firmly  welded  to  the  steel  after  the  explosion. 

Figure  6  is  a  block  diagram  of  the  apparatus.  The  Mk  88  electronic 
firing  switch  supplied  two  simultaneous  pulses,  each  of  about  1000  v, 
from  a  1-pf  capacitor.  One  pulse  fired  the  EBW  detonator;  the  other 
drove  the  primary  of  an  air  core  transformer.  This  transformer,  with 
about  a  20-to-l -turns  ratio,  stepped  up  the  voltage  of  the  pulse  from  the 
Mk  88  firing  switch  sufficiently  to  trigger  the  spark  gap  switch;  thus  the 
EBW  detonator  and  the  spark  gap  switch  were  fired  simultaneously. 

The  tantalum  strip  load  formed  the  resistance  R2  of  one  arm  of  the 
bridge  shown  in  Figure  7.  The  resistive  voltage  drop  across  the  tanta¬ 
lum  strip  is  approximately  equal  to  the  inductive  voltage  drop  produced 
by  L  across  the  same  strip.  A  section  of  the  upper  copper  plate, 
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FIGURE  6.  Block  Diagram  of  Magnetic  Field  Comp 


sufficiently  long  to  produce  an  inductance  equal  to  that  of  the  inductance 
of  the  tantalum  strip,  formed  another  arm  of  the  bridge  R1  and  L.  The 
other  parts  of  the  bridge  were  resistors  (see  Figure  5)  shown  schemati¬ 
cally  in  Figure  7.  The  output  of  the  bridge  was  conducted  by  a  twisted 
pair  of  wires  to  an  attenuator  and  pulse  transformer  placed  about  6  ft 
from  the  field  compressor  and  outside  the  firing  chamber.  The  pulse 
transformer  was  necessary  to  prevent  the  large  ground  currents  flowing 
in  the  flux  compressor  from  upsetting  the  operation  of  the  oscilloscopes 
used  for  recording  the  voltage  developed  in  the  bridge  of  the  field  com¬ 
pressor. 

The  pulse  transformer  was  constructed  from  two  ferrite  cores,  each 
in  the  form  of  a  letter  E.  Copper  wire  was  wound  on  the  center  arm  of 
the  E's  and  the  two  sections  were  then  clamped  together  with  electric 
tape.  The  transformer  had  a  1-to- 1-turns  ratio.  A  test  of  the  trans¬ 
former  showed  that  it  would  pass  a  100-fisec  pulse  without  significant 
distortion  as  long  as  the  input  pulse  was  not  more  than  5  v.  A  voltage 
divider  was  therefore  connected  to  the  input  of  the  pulse  transformer  to 
keep  the  voltage  below  this  level. 

Oscilloscopes  and  a  high  voltage  power  supply  for  charging  the  capa¬ 
citor  bank  were  placed  in  a  firing  shack  and  connected  to  the  capacitor 
bank  and  field  compressor  by  100  ft  of  high  voltage  coaxial  cable. 

To  test  the  apparatus,  the  capacitor  bank  was  first  charged  to  3000  v 
and  then,  without  a  detonator,  the  firing  circuit  was  armed  and  fired, 
thus  triggering  the  spark  gap  switch  and  the  oscilloscopes.  The  result 
is  shown  in  Figure  8,  which  is  a  photograph  of  an  oscilloscope  trace. 

A  detonator  was  placed  on  the  field  compressor  and  the  capacitor  bank 
was  again  charged  to  3000  v.  This  time  when  the  firing  switch  was 
closed,  the  spark  gap  switch  discharged  the  capacitors  into  the  field 
compressor  and  the  PETN  sheet  explosive  forced  the  two  copper  plates 
together.  The  oscilloscope  trace  of  this  test  is  shown  in  Figure  9.  The 
maximum  current  was  about  2.8  times  that  obtained  without  the  explosive 
charge.  Current  increase  stopped  about  34  y.sec  after  field  compression 
began.  The  detonation  wave  in  the  PETN  sheet  travels  about  7  mm/Vsec, 
and  field  compression  stopped  after  the  detonation  wave  had  advanced 
about  9.4  in.  Why  the  field  compression  stopped  then  is  not  known. 

Figure  10  is  a  photograph  taken  of  the  compressor  after  the  test. 

The  figure  shows  that  the  upper  and  lower  copper  strips  were  still  con¬ 
tinuous,  although  there  was  a  tear  in  the  upper  strip.  The  upper  strip 
had  become  welded  to  the  lower  strip  along  the  edges,  and  it  also 
exhibited  bulges  as  if  some  fluid  had  been  trapped  under  it.  Calculations 
indicate  that  the  magnetic  field  pressure  was  only  a  few  hundred  pounds 
per  centimeter  of  compressor  length. 
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FIGURE  9.  Field  Compressor  Output  With 
Explosive.  (Time  scale:  20  ^sec/crn) 


After  9.4  in.  of  compression,  not  enough  energy  would  have  been 
transferred  to  the  tantalum  strip  to  heat  it  sufficiently  to  increase  its 
resistance  significantly,  and  so  (5  of  equation  (8)  must  have  remained 
near  its  initial  value  of  0.1.  From  equation  (8),  with  X  =  9.4,  Xq  =  15, 
and  (3  =  0.1 


which  is  close  to  the  measured  value  of  2.8.  The  size  of  the  tantalum 
strip  was  based  on  a  compression  ratio  of  20.  If  further  experiment  s 
show  that  a  compression  ratio  this  large  cannot  be  obtained,  the  resist¬ 
ance  of  the  tantalum  strip  load  should  be  increased. 
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SEMICONDUCTOR  PHYSICS 


ACCURATE  MEASUREMENT  OF  SMALL 
ROTARY  DISPERSION 

by 

C.  J.  Gabriel  and  H.  Piller 
Infrared  Division 


INTRODUCTION 

The  study  of  interband  Faraday  rotation  in  semiconductors  requires 
the  transmission  of  electromagnetic  energy  through  highly  absorbing 
materials.  In  order  to  have  sufficient  transmitted  energy  to  make  the 
measurements,  often  samples  must  be  used  having  thicknesses  of  the 
order  of  l/q,  where  n  xs  the  absorption  coefficient.  The  rotations 
induced  by  such  thm  samples  can  be  less  than  one  minute.  This  report 
discusses  how'  the  addition  of  phase- sensitive  detection  to  an  appa¬ 
ratus  described  in  a  previous  report*  allows  the  measurement  of 
small  rotations. 

MODIFICATION  OF  THE  APPARATUS 

Priefly,  the  apparatus  consists  of  a  high-pressure  zenon  arc,  focusing 
optics,  a  510-cps  light  chopper,  monochromator,  Gian- Thompson 
polarizer,  sample  cell.  Wollaston  prism  beam  splitter,  two  photodetec¬ 
tors.  a  bridge  circuit,  ana  suitable  amplifiers.  Details  of  the  apparatus 
were  given  previously.^ 

The  modification  was  made  by  replacing  the  Hewlett-Packard  Model 
415A  amplifier  by  a  Princeton  Applied  Research  Model  JB-5  lock-in 
amplifier.  Because  the  gain  of  the  PAR  JB-5  is  somewhat  less  than 
the  gain  of  the  HP  415A,  a  Tektronix  122  preamplifier  was  used.  Since 
the  Tektronix  122  is  double-ended,  the  coupling  transformer  was  elimi¬ 
nated  and  capacitive  coupling  was  used  instead.  A  parallel  LC  circuit 
turvd  to  the  chopping  frequency  was  placed  at  the  input  of  the  preampli¬ 
fier;  its  purpose  was  to  decrease  the  loading  effect  on  the  tur.  d  stages 


*C.  J.  Gabriel,  V.  A.  Patton,  and  H.  Piller,  NAVWEPS  Report  8226, 
Quarterly  Report:  Foundational  Research  Projects.  October- December 
1964,  p,  69. 

2Ibid. 
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of  the  lock-in  amplifier  by  the  pickup  of  60-cps  voltage  from  the  line. 

The  reference  for  this  amplifier  was  obtained  by  mounting  a  flashlight 
lamp  and  solar  cell  on  opposite  sides  of  the  chopper  disc.  The  voltage 
obtained  from  the  solar  cell  was  sufficient  to  drive  the  reference  channel 
of  the  lock-in  amplifier. 

Three  advantages  result  from  the  use  of  the  phase- sensitive  amplifier: 

1.  An  improvement  is  obtained  in  signal-to-noise  ratio  because  of 
the  narrow  bandpass  of  the  lock-in  amplifier.  This  effect,  being  well 
known, 3  will  not  be  discussed  here. 

2.  A  bipolar  signal  is  obtained.  The  sign  of  the  signal  indicates 
which  detector  is  more  stongly  irradiated,  thus  making  the  balancing  of 
the  bridge  easier  and  indicating  the  sign  of  the  rotation. 

3.  A  more  linear  response  is  obtained.  This  allows  the  accurate 
measurement  of  small  rotation,  regardless  of  phase  imbalance  of  the 
bridge  network.  This  point  will  be  discussed  in  detail. 


THEORY  OF  MEASUREMENT 


Assume  that  the  signal  voltages  from  the  photodetector  are  out  of 
phase  by  the  angle  4>  and  that  4>  is  independent  of  the  signal  amplitudes 
vj  and  V£.  Let  the  phase  4  of  the  reference  signal  be  measured  from 
the  bisector  of  the  phase  difference  4>.  (See  Figure  1.)  The  components 
of  the  difference  signal  vs  and  the  sum  signal  V  that  are  in  phase  with 
the  reference  are  then  given,  respectively,  by  the  expressions 


v  (a  +  0)  =  (v 
s  l 


d>  <j> 

v^'cosfccos-j  -  (Vj  +  v^)sin$sin  — 


V(a  +  0)  =  (v^  +  v^)cos  $cos  j  -  (v^ 


v^)  sin  $  sin  ~ 


where  the  signal  amplitudes  are  given  by  the  expressions 

Vj  =  IKjffjkj  cos2  (a  +  0)  (3) 

v ^  s  IK^cr^k^  sin2  (a  +  0)  (4) 

^R.  H.  Dicke,  Rev.  Sci.  Instr.,  Vol.  17,  p.  268  (1947). 
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FIGURE  1.  Relative  Phases  of  the  Signal  Voltages  vj,  V2» 
and  the  Reference  Channel 


in  which  I  is  the  source  intensity,  Kj  and  K2  are  factors  depending  on 
the  bridge  parameters,  <rj  and  <T2  are  factors  depending  on  the  detector 
responses,  kj  and  k2  are  factors  to  take  into  account  differences  in 
absorption  in  the  two  beams,  a  is  the  angle  between  the  planes  of  polari 
zation  of  the  polarizer  and  analyzer,  and  0  is  the  rotation  of  the  plane 
of  polarization  caused  by  rotary  dispersion  in  the  sample. 


The  expressions  (1)  and  (2)  for  vB  and  V  should  be  compared  with  the 
corresponding  expressions  for  the  amplitudes  of  v8  and  V  that  are  ob¬ 
tained  if  phase- sensitive  detection  is  not  used: 


v 

s 


V 


+  vz> 


Zvivi 


(v,  -  v2> 


<1 


2v1v2 


(Vj  +  v2) 


<1 


cos 


cos 


♦)1 

♦)1 


(5) 


(6) 


It  is  evident  that  equation  (1)  is  linear  in  0  in  the  region  where  v^  « 
whereas  equation  (5)  is  not.  If  V  were  obtained  by  measuring  vj  and  V2 
separately,  by  covering  eacn  detector,  then  equation  (6)  would  become 
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V  =  +  V£.  It  is  assumed  that  V  is  determined  from  a  single  meas¬ 

urement. 

For  values  of  0  less  than  5  degrees,  an  expansion  to  second  order 
in  0  can  be  used. 

6  T  An  <t>  / 

vg(a  +  0)  «  -incostcos  ^  I  sin  2(3  -  +  tan  $  tan  —  ll 

+  2  ^1  -  ~—tan  $  tan  ^0  cos  2(3 

,  TT 

where  a  =  —  +  B 
4 

2n  =  KjO-jkj  +  K2a2k2 

2An  =  KjtTjkj  - 

A  similar  expression  is  obtained  for  V(a  +  0)  except  that  An  and  n  are 
interchanged.  The  system  can  generally  be  adjusted  so  that  |3,  An/n, 
and  4>  «  1.  Equation  (7)  and  the  similar  expression  for  V(a  +  0)  reduce 
to: 

vg(a  +  0)  «  -incos  |^20  +  20  -  +  *  tan  (8) 

V(q  +  0)  «  -incos  ijip—  -  1  +  #tan|(~  -  2pj 

+  2e(f 

where  terms  to  third  order  in  small  quantities  have  been  dropped  and 
it  has  been  assumed  that  $  <  ir/2. 


The  Faraday  rotation  is  obtained  from  the  measured  values  of  vg  and 
V  by  either  of  the  following  combinations:  Ri  =  (v£  -  v~ )/(V+  +  V")  or 
R.2  =  (v*/2V+)  -  (vg/2V"),  where  the  superscripts  (+)  and  (-)  indicate 
measurements  made  with  forward  and  reversed  magnetic  fields,  respec¬ 
tively.  From  equations  (8)  and  (9),  by  dropping  small  quantities  of 
third  order, 
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(11) 


where  20  =  0+  -  0 
2  AH  =  H+  +  H" 
2H  =  H+  -  H“ 
21  =  I+  +  I" 

V  V  V 

2AI  =  I+  -  I' 
v  v  v 


with  similar  definitions  for  Iy  and  Aly,  and 


A 


B 


The  subscripts  v  and  V  designate,  respectively,  the  average  value  of  the 
source  intensity  when  either  v8  or  V  is  being  measured,  and  H+  and  H” 
are  the  average  values  of  the  forward  and  reversed  magnetic  fields, 
respectively. 


DISCUSSION 

Generally  0  is  obtained  from  equation  (10)  or  equation  (11)  by  neglect¬ 
ing  all  unmeasured  quantities;  thus,  6  =  K\/Z  or  0  s  R2/2.  In  addition 
to  this  error  of  negligence,  an  error  is  introduced  because  of  the  error 
involved  in  measuring  v8  and  V.  For  the  moment  this  error  will  be 
neglected.  The  value  of  0  obtained  in  this  way  is  the  sum  of  the  rotation 
08  induced  by  transmission  through  the  sample  and  the  rotation  0m 
caused  by  reflection  from  mirrors  located  in  the  magnetic  field.  The 
mirrors  are  located  in  the  magnetic  field  because  the  magnet  used  at 
NOLC  does  not  have  axial  access.  The  value  of  0m  must  be  obtained 
from  a  separate  measurement.  For  simplicity,  consider  equation  (10); 
the  treatment  of  equation  (11)  is  similar.  Equation  (10)  has  the  form 

R  =  26(1  +  p)  +  q  (12) 
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A  similar  expression  holds  for  the  measurement 


where  9  =  9a  +  9  . 

s  m 

of  9m;  thus 


R 

m 


29  (1  +  p  )  +  q 
m  m  m 


(13) 


From  equations  (12)  and  (13),  an  expression  for  the  error  in  the  meas¬ 
ured  value  of  9g  can  be  obtained. 

9,  -  (f  ■  -q  ■  -P0,  +  <Pm  -  P>9m  +  ^  ’  «>  <‘4> 

Because  of  the  random  nature  of  the  quantities  AIv/lv  and  AH/H,  the 
factors  pm  -  p  and  qm  -  q  are  not  necessarily  zero. 

The  first  term  of  equation  (14)  is  of  the  order  of  (AH/H)  (AIv/Iv)98, 
and  the  second  term  is  of  the  order  of  (AIv/lv)  (AH/H)9m.  The  magnitude 
of  the  third  term  depends  on  the  method  used  to  balance  the  system.  In 
principle,  if  the  system  is  adjusted  until  v8(a)  ~  0,  then  q  =  qm  =  0  (see 
equation  8).  The  difficulty  with  this  adjustment  is  that  the  magnetic 
field  can  cause  a  movement  of  the  optical  components  and  thus  produce 
a  displacement  of  the  light  beams  with  respect  to  the  detectors,  resulting 
in  an  imbalance,  (A ft/tt)^,  of  the  system.  Since  the  imbalance  does  not 
depend  on  the  direction  of  the  magnetic  field,  it  does  not  appear  as  a 
rotation,  but  only  contributes  to  qm  -  q.  This  contribution  is  of  the 
order  of  (An/n)^  ( AIv/ly).  Alternatively,  the  system  can  be  adjusted 
until  the  condition  v^  +  v  ~  -  0  holds.  Under  this  condition 
qm  -  q  (9m  +  9s)(Alv/lv  +  AH/H)  Aly/ly.  The  second  balancing  method 
generally  gives  the  best  results. 

Instead  of  making  fine  adjustments  in  the  system,  an  alternative  pro¬ 
cedure  is  available.  If  the  electronics  of  the  detection  system  can  be 
arranged  so  that  the  measurements  of  vs  and  V  for  a  given  field  direc¬ 
tion  take  place  simultaneously,  or  nearly  so,  then  equation  (11)  can  be 
used  virtually  without  error  introduced  by  neglecting  unmeasured  quanti¬ 
ties.  This  comes  about  because  under  this  condition  Iv  «  Iy  and 
AIV  «  Aly  so  that  B  *  0  and  A  «  1 . 

In  contrast  to  phase  detection,  fine  adjustement  of  the  system  is  not 
possible  with  amplitude  detection  because,  as  can  be  seen  from  equation 
(5),  v8  will  be  linear  only  if  the  system  is  unbalanced — for  instance,  by 
making  0  >  $.  This  gives  an  error  in  9  proportional  to  $AIv/lv.  How¬ 
ever,  the  alternative  procedure  mentioned  above  is  still  available. 

If  the  system  is  such  that  AH/H  w  0.01,  AIv/lv  »  0.1,  and  9g  0m, 
then  the  error  of  negligence  will  be  of  the  order  of  9m  X  10"^.  The 
polar  Faraday  rotation  on  the  mirrors  in  a  field  of  20  kG  is  typically 
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less  than  0.001  radian.^  Thus  the  error  in  0  will  be  of  the  order  of  10"^ 
radian.  This  error  is  of  the  same  order  as  that  introduced  by  errors 
in  the  measurement  of  vg  and  V;  that  is,  under  optimum  conditions 
Avg/V  «  10“6,  where  Avg(  the  rms  noise  voltage  in  vg,  is  typically  of 
the  order  of  a  microvolt  and  V  is  typically  of  the  order  of  a  volt.  The 
values  of  Avg  and  V  quoted  above  are  obtained  by  using  PbS  detectors, 
time  constants  of  several  seconds,  and  samples  thin  enough  so  that  the 
transmitted  intensity  is  about  0.3  of  the  incident  intensity. 

A  severe  restriction  on  this  theoretical  sensitivity  of  10"^  to  10"^ 
radian  arises  if  the  apparatus  is  not  mechanically  stable.  If  the  appara¬ 
tus  is  rigidly  constructed,  with  attention  being  paid  to  the  isolation  of 
the  system  from  the  surroundings,  then  movements  of  the  optical  com¬ 
ponents  can  be  negligible.  However,  instabilities  in  the  intensity  distri¬ 
bution  of  the  arc  itself  will  still  be  present.  Since  such  instabilities 
introduce  an  error  into  6  of  the  order  of  6(Afi)/ft,  where  5(Afl)  is  the 
change  in  Afi  during  the  time  interval  between  the  measurements  of  v  * 
and  v  ",  the  sensitivity  can  be  reduced  by  several  orders  of  magnitude. 
For  this  reason,  the  detectors  should  be  positioned  so  that  movement 
of  the  arc  affects  both  detectors  equally.  Moreover,  it  is  advantageous 
to  sacrifice  the  signal-to-noise  ratio  of  the  detectors  and  to  use  detec¬ 
tors  which  are  large  enough  for  the  beams  to  remain  within  the  sensitive 
area  for  small  movements  of  the  arc. 

SUMMARY 

It  has  been  shown  that  the  utilization  of  phase-sensitive  detection  with 
a  double-beam  null  technique  provides  an  increase  in  the  linearity  of 
the  response  of  the  system;  moreover,  that  this  increase  in  linearity 
allows  the  measurement  of  small  values  of  Faraday  rotation  independent 
of  phase  imbalance  of  the  bridge  circuit,  residual  magnetization  of  the 
electromagnet,  and  source  intensity  fluctuations,  provided  these  quanti¬ 
ties  are  reasonably  small.  However,  instabilities  in  the  intensity  dis¬ 
tribution  of  the  arc  limit  these  advantages  to  such  an  extent  that  the 
detector  signal-to-noise  ratio  should  be  sacrificed  by  using  detectors 
large  enough  so  that  the  beams  remain  within  the  sensitive  areas  for 
small  movements  of  the  arc.  If,  by  the  use  of  detectors  other  than  PbS, 
the  ratio  Avg/V  can  be  reduced  to  below  10’^,  then  the  apparatus 
should  be  modified  in  order  to  fully  utilize  the  detectors.  The  system 
should  be  so  designed  that  no  mirrors  are  located  in  the  magnetic  field 
and  the  ratio  vg/V  is  measured  directly. 
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ELLIPTICITY  IN  MULTIPLY  REFLECTED 
FARADAY  ROTATION 

by 

C.  J.  Gabriel  and  H.  Piller 
Infrared  Division 


INTRODUCTION 

The  presence  of  multiple  internal  reflections  can  cause  difficulty  in 
the  accurate  measurement  of  Faraday  rotation  and  ellipticity  in  slab¬ 
shaped  samples  of  materials  for  which  the  product  of  the  absorption 
coefficient  tj  and  the  thickness  d  satisfies  the  condition  Tjd  &  1.  In  a 
previous  report,^  a  model  was  discussed  whereby  multiply  reflected 
Faraday  rotation  could  be  conveniently  treated.  In  this  report  the  model 
will  be  extended  to  include  the  effects  of  internal  polar- reflected  Faraday 
rotation.  The  primary  purpose  of  the  report,  however,  is  to  provide 
simple  relationships  between  the  observed  rotation  and  ellipticity  and 
the  rotation  produced  by  a  single  passage  of  the  light  beam  through  the 
sample.  Two  extreme  cases  will  be  considered:  (1)  that  of  precise 
coherence  among  the  multiply  reflected  beams,  and  (2)  that  of  complete 
incoherence  among  the  multiply  reflected  beams.  The  first  case  is 
especially  useful  in  the  treatment  of  experiments  that  use  a  laser  as  a 
light  source  and  thin,  highly  polished  samples.  The  second  case  applies 
in  low- resolution  spectroscopy  and  in  experiments  utilizing  samples  that 
are  rough  compared  to  the  wavelength  of  the  light  used. 

THE  MODEL 

The  following  model  is  used  to  discuss  the  effects  of  multiple  reflec¬ 
tions. 

1.  A  collimated,  linearly  polarized  light  beam  is  normally  incident 
upon  a  plane  parallel  slab. 

2.  The  amplitude  of  the  transmitted  beam  of  light  is  the  sum  of  the 
amplitudes  of  the  light  beams  produced  by  partial  internal  reflec¬ 
tion  on  the  sample  surfaces. 

3.  The  effect  of  the  magnetic  field  is  introduced  phenomenologically 
as  a  rotation  of  the  plane  of  polarization  upon  each  passage  through 


^H.  Piller  and  C.  J.  Gabriel,  NAVWEPS  Report  8812,  Quarterly 
Report:  Foundational  Research  Projects,  January-March  1965,  p.  27. 
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the  sample  and  upon  each  internal  reflection  from  the  surfaces. 
The  ellipticity  produced  in  a  single  passage  through  the  sample  is 
neglected. 


Only  a  brief  treatment  of  this  model  will  be  given  here. 


After  passing  through  an  analyzer  rotated  by  an  angle  a  with  respect 
to  the  polarizer,  the  amplitude  of  the  q  times  reflected  beam  is  given  by 


=  t1t2(-r1r2)q  e  ^2q+1^cos  [(2q  +  1)0  +  2q$  +  <$>'  -  a] 


(1) 


where  tj,  t2,  rj,  and  r2  are  the  Fresnel  coefficients  for  normal  inci¬ 
dence: 


1  -  n  +  ik 
1  1  +  n  -  ik 


r.  = 


t.  = 


1  n  +  1  -  ik 


r„  = 


t,  - 


n  -  1  -  ik 
1  +  n  -  ik 

2(n  -  ik) 
n  +  1  -  ik 


(2) 


in  which  n  and  k  are  the  real  and  imaginary  parts  of  the  complex  index 
of  refraction  of  the  sample,  and  n  =  1,  k  =  0  in  the  surrounding  medium. 
The  complex  phase  shift  6  is  given  by  6  =  2w(n  -  ik)d/X,  where  X  is  the 
wavelength  in  vacuum.  The  factor  exp[-i(2q  +  1)6]  in  equation  (1)  takes 
into  account  the  phase  shift  and  absorption  occurring  in  the  bulk  material. 
The  angle  0  is  the  Faraday  rotation  for  a  single  passage;  $  is  the  polar 
Faraday  rotation  that  takes  place  upon  each  internal  reflection;  and 
is  any  magneto  rotation  that  might  take  place  on  passage  through  both 
surfaces.  Although  is  expected  to  be  zero,  it  is  included  for  com¬ 
pleteness. 

The  intensity  of  the  emerging  beam  will,  in  general,  be  a  periodic 
function  of  the  analyzer  position  a.  The  value  of  a  for  which  the  intensity 
is  a  relative  maxin  um  is  the  observed  rotation  and  will  be  designated 
by  v*  The  ratio  of  the  minimum  intensity  to  the  maximum  intensity  gives 
the  ellipticity  in  the  case  of  complete  coherence,  and  a  measure  of  the 
depolarization  in  the  case  of  complete  incoherence. 


COMPLETELY  COHERENT  BEAMS 

The  transmitted  amplitude  for  coherent  beams  is  given  by  the  summa¬ 
tion  of  the  amplitudes  obtained  from  equation  (1).  Hence, 


with  d  the  sample  thickness  and  X  the  wavelength.  The  quantities  a  and 
b  are  the  real  and  imaginary  parts,  respectively,  of  rjr;>;  R  is  the 
reflectivity;  and  r\  is  the  absorption  coefficient.  The  parameter  K  (see 
Appendix  A)  is  independent  of  a  and  thus  an  explicit  form  is  not  necessary 
here. 


The  condition  d(rr*)/da  *  0  yields  the  expression 


tan  Zr  * 


(1  -  c^)  sin  2© 

(1  +  c^)  cos  2©  7  c  j 


where  r  *  y  +  +  -  +'.  This  expression  relates  the  observed  rotation  y 
to  the  Faraday  rotation  9  and  the  polar  reflection  Faraday  rotation  ♦, 
since  can  be  expected  to  be  aero.  The  partial  inversion  of  equation 
(10)  yields 
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* 


sin  2(6  -  4)  = 


c  tan  2r 


[(1  -  c  )2  +  (1  +  c  )2tan22r]i 


where  tan  24  =  ((1  +  c2)/(l  -  c2)]tan2I\ 

The  transmitted  beam  becomes  eliiptically  polarized  upon  passing 
through  the  sample.  The  ellipticity  arising  in  this  model  is,  of  course, 
not  related  to  the  intrinsic  ellipticity,  which  has  been  neglected  here, 
but  is  caused  by  the  multiple  reflections.  The  ellipticity  C  is  defined  by 
€2  -  ^min/^max'  w^ere  and  are,  respectively,  the  minimum 

and  maximum  intensities  obtained  from  equation  (4)  as  a  is  varied  through 
180  degrees.  Since  occurs  when  A  =  T,  then  I^j,  occurs  when 

A  =  T  +  ir/2.  Hence,  from  equation  (4) 

cos  2(6  -  T)  +  c  j  cos  2T  +  c^  cos  2(6  +  T) 

p  _______  (12j 

where  p  =  (1  -  €2)/(l  +  €2).  In  principle,  equations  (12)  and  (10)  can  be 
used  to  obtain  8  and  6  from  measured  values  of  c2  and  y  if  the  condition 

*  0  holds. 

By  treating  sin  26  and  cos  26  as  independent  variables,  equations  (10) 
and  (12)  can  be  solved  simultaneously,  to  yield 


sin  26  s 


(1  tca)2-c|2 

1  -  C, 


p  sin  2T 

1  4  c^  -  pc  j  cos  2T 


cos  26  - 


p(l  +  c^) cos  2T  -  c  j 
1  +  c^  -  pc  j  cos  2T 


From  the  identity  sin2  26  4  cos2  26  s  1,  it  follows  directly  that 

4c  .  c*  Ti 
^2  Cl  2 

p  =  If  — - - j  sin  2T  ( 1 

o  -  c2r 

*  m 

Thus,  r  can  be  determined  from  the  measured  value  of  p.  Then,  using 
equation  (10),  6  can  be  obtained.  If  the  conditions  ♦  *  0  and  4'  »  0  hold, 
then  equation  (15)  can  be  used  to  check  the  validity  of  the  model  by  re¬ 
lating  the  measured  quantities  y  and  p. 
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COMPLETELY  INCOHERENT  BEAMS 


In  the  case  of  incoherent  beams,  the  intensity  of  the  transmitted  beam 
is  given  by 


=  2  T  T  * 


,T0  <>  * 


Evaluation  of  equation  (16)  yields 


tt*  =  K'{l  +  c  2  -  2c2cos40  +  (1  -  c ^ )  [c o s  2(6  -  A) 


c^  cos 


2(0  +  A)]l 


where  K'  (see  Appendix  A)  is  independent  of  a.  The  condition 
d(TT*)/da  =  0  yields 


1  -  c. 


tan  20  =  - — - tan  2T 

1  +c2 


This  equation  relates  0  and  <j>  to  \  since  the  condition  <f>*  —  0  is  expected 
to  hold. 


The  transmitted  beam  is  not  elliptically  polarized  but,  rather,  par¬ 
tially  depolarized.  The  polarization  factor  (Imax  -  ^min^/t^max  +  *min) 
has  the  same  form  as  the  quantity  p  used  above  to  describe  the  ellipticity. 
The  same  symbol,  p,  will  be  used  for  the  polarization  factor.  From 
equation  (17) 


1  -  c. 


1  +  c^  -  2c2  cos  40 


[cos  2(0  -  T)  -  c  cos  2(0  +  T)] 


The  insertion  of  equation  (18)  into  equation  (19)  yields 


4c 

(1  + 


sin22T  2 


Equation  (20)  gives  F  in  terms  of  the  measured  quantity  p.  This  value 
of  r,  when  substituted  in  equation  (18),  yields  a  value  of  0.  Thus,  from 
the  measured  values  of  p  and  y,  both  9  and  4>  can  be  obtained  if  the  con- 


dition  of  =*  0  holds.  Under  the  conditions  4>  “  0  and  =*  0,  the 
applicability  of  the  model  can  be  checked  by  using  equation  (20)  to  relate 
the  measured  values  of  y  and  p. 


DISCUSSION 

Although  the  rotation  8  can  be  determined  from  a  measurement  of  y, 
assuming  4»  =  d>*  —  0,  photometric  techniques  generally  do  not  measure 
y  independently  of  p,  so  either  a  measurement  or  calculation  of  p  is 
necessary.  For  example,  in  a  technique  used  at  this  Laboratory,^  the 
transmitted  light  beam  is  split  into  two  perpendicularly  polarized  com¬ 
ponents.  When  the  system  is  properly  balanced,  the  quantity  that  is 
measured  is  (I i  -  I2)/ (^1  +  l2)>  w^ere  Ij  and  I2  are  the  intensities  of  the 
two  components  of  the  light  beam.  Writing 

1  =1  ,cos2(a  -  y)  +  I  sin2  (a  -  y)  (21) 

i  it.i  t  min 

J2  =  Lv8i^a  -  v)  +  IrY,incos2(a  -  Y)  (22) 

4  max  min 

gives 


f-;  f  =  p  cos  2(a  -  y)  (23) 

J1  +I2 

Generally,  the  system  is  adjusted  so  that  a  =  t/i.  Under  this  condition 
Ul  -  I2)/(Ii  +  l2)  =  p  sin  Zy.  A  separate  measurement  of  p  can  be  made 
if  the  system  is  adjusted  until  a  =  y,  that  is,  for  maximum  Ij  -  1^1  then 
(If  -  1 2 )/ (I l  +  I2)  =  P»  However,  if  p  —  1,  a  measurement  of  p  made  in 
this  way  can  be  more  inaccurate  than  simply  to  set  p  =  1  initially  and 
then  make  successive  approximations  until  consistency  is  reached  by 
using  the  relationships  among  p,  y,  and  8  and  the  measured  value  of 
psin2y.  For  small  values  of  8,  equations  (15)  and  (20)  indicate  that 
p  =  1,  and  y  is  obtained  directly  from  the  measurement. 

In  order  to  use  the  expressions  for  the  rotation  and  ellipticity,  it  is 
necessary  to  know  the  values  of  cj  and  C£.  This  requires  the  accurate 
knowledge  of  n,  k,  d,  and  \,  particularly  if  interference  effects  are 
observed.  For  example,  suppose  that  ®  is  only  known  to  within  one  per¬ 
cent;  then  the  error  in  will  be  0,0 l(4irnd/\)  =*  O.lnd/X,.  Unless 

C.  J.  Gabriel,  V.  A.  Patton,  and  H.  Piller,  NAVWEPS  Report  8226, 
Quarterly  Report:  Foundational  Research  Projects,  October-December 
1964,  p.  69. 
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V  is  large,  nd/X  can  be  greater  than  one,  so  the  condition  A®  «  tt  does 
not  hold,  and  cj  cannot  be  obtained  accurately.  When  interference  effects 
are  not  observed,  the  situation  is  improved  since  cj  does  not  appear  in 
the  expressions  describing  this  case  and  C£  does  not  involve  the  phase  $. 

Even  if  Cj  and  eg  are  not  known  accurately,  the  expressions  derived 
are  still  useful  for  the  interpretation  of  experiments  in  which  the  rota¬ 
tion  or  the  ellipticity  is  measured  as  a  function  of  d,  X,  or  the  magnetic 
field.  Moreover,  by  using  these  expressions,  independent  measurements 
of  both  v  and  p  provide  a  check  on  the  consistency  of  the  values  of  cj 
and  c 2*  and  thus  a  check  on  the  accuracy  of  the  correction  necessary  to 
obtain  the  intrinsic  rotation. 

APPENDIX  A— VALUES  OF  K  AND  KJ 

In  order  that  equations  (4)  and  (17)  can  be  conveniently  used  for  the 
interpretation  of  experiments  involving  the  magneto-absorption  of 
linearly  polarized  light,  the  explicit  forms  of  K  and  K'  are  given  here. 

i  y  y  i  1 

K  =  Te  ^[[l  +  2c  j(  J  +  c^)  coc  23  +  c  ^  +  2c ^  cos  40  +  c^  ]  (24) 

K'  =  Te”1ld[2(l  -  c2)(l  -  2c2cos40  +  c2)]  1  (25) 

where  T  =  1 1 j t^  |  From  equation  (2),  Im(tjt2)  =  b  and 


OPTICAL  CONSTANTS  OF  GERMANIUM  IN  SPECTRAL 
REGIONS  FROM  0.5  TO  3.0  eV 
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INTRODUCTION 

The  interaction  of  electromagnetic  radiation  with  materials  has  long 
been  recognized  as  a  means  of  acquiring  information  on  the  basic  proper¬ 
ties  of  the  materials.  For  those  instances  where  the  Maxwell  theory  of 
electromagnetism  can  be  expected  to  hold,  the  measurements  are  usually 
presented  in  the  form  of  optical  (or  dielectric)  constants.  The  word  con¬ 
stants  would  better  be  replaced  with  the  term  coefficients,  inasmuch  as 
the  constants  vary  considerably  as  a  function  of  wavelength,  material 
source,  and  techniques  of  measurement.  It  is  this  variation,  of  course, 
which  provides  the  spectroscopic  information  that  is  sought.  The  meas¬ 
urement  of  the  optical  properties  of  opaque  but  specular  materials  has 
had  increased  activity  in  recent  years.  The  incentive  has  come  mainly 
from  the  need  for  more  experimental  details  on  the  band  structure 
parameters  of  semiconductors.^  Information  concerning  the  band  struc¬ 
tures  of  the  metals  has  also  been  gained  from  studies  of  the  optical 
properties.® 

Three  principal  techniques  have  been  used  for  measuring  the  complex 
dielectric  constants  €  =  €j  +  iCg  in  the  near  infrared,  visible,  and  ultra¬ 
violet  regions.  The  technique  that  has  provided  the  most  data  in  recent 
years  has  been  the  measurement  of  reflectance  (Rq)  with  the  radiation 
incident  nearly  normal  to  the  sample  surface.  Much  of  this  work  has 
been  done  without  an  absolute  reflectance  measurement  but  in  arbitrary 
units.  Interpretation  of  such  spectra  depends  upon  additional  information 
as  to  whether  the  reflectance  is  dominated  by  the  real  or  the  imaginary 
part  of  the  dielectric  constant. 

Using  this  technique,  Tauc^  and  Cardona^  have  made  considerable 
contributions,  providing  high  resolution  data  over  a  wide  spectral  region. 

^J.  C.  Phillips,  Phys.  Rev.,  Vol.  133,  p.  A  452  (1964). 
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For  a  recent  review,  see  M.  Suffczynski,  Phys.  Stat.  Sol.,  Vol.  4, 

p.  3  (1964). 

a 

For  a  review,  see  J.  Tauc,  Proceedings  of  the  International  Confer¬ 
ence  on  the  Physics  of  Semiconductors,  Exeter  (1962),  p.  333;  also 
S.  Lukei  and  A.  Schmidt,  op.  cit.,  p.  389. 

*^M.  Cardona,  J.  Appl.  Phys.,  Vol.  32,  p.  958  (1961). 


Much  of  their  work  was  concerned  with  spectral  regions  where  €£*  the 
imaginary  part  of  the  dielectric  constants,  is  dominant  and  thus  Rq  has 
spectral  features  similar  to  €£•  Consequently,  from  Rq  one  can  gain 
an  indication  of  the  relative  shape  of  the  density-of- states  function. 

This  technique  is  not  so  suitable,  however,  in  the  regions  where  the 
dominance  of  either  or  € 2  is  indeterminant. 

The  measurement  of  the  absolute  reflectance  at  normal  incidence  has 
been  carried  out  by  several  workers.  Philipp  and  Taft**  applied  the 
so-called  Kramers-Kronig* ^  analysis,  which  is  applicable  if  sufficient 
spectral  range  has  been  covered.  Given  the  reflectance  over  a  sufficient 
spectral  region,  the  real  and  imaginary  parts  of  the  dielectric  constant 
can  be  calculated  for  that  region.  The  principal  requirement  is  to  insure 
that  all  of  the  necessary  information  has  been  gathered. 

The  technique  of  ellipsometry,  which  also  gives  the  real  and  imaginary 
parts  of  the  dielectric  constant,  has  been  used  to  study  germanium. 
Archer  applied  this  technique  to  examine  the  optical  properties  of  ger¬ 
manium  in  the  visible  region.*-^  His  results  are  in  good  agreement  with 
those  of  Philipp  and  Taft  and  with  more  recent  measurements  by  Philipp 
and  Ehrenreich*^  which  were  obtained  from  the  K-K  approach.  Archer's 
results  are  also  in  reasonable  agreement  with  recent  measurements  by 
Marple  and  Ehrenreich,*®  who  also  applied  ellipsometric  techniques, 
although  Archer's  do  not  resolve  certain  spectral  features. 

Photometric  techniques  have  also  been  used  for  determining  the  opti¬ 
cal  constants.  Avery*®  and,  more  recently,  Lindquist  and  Ewald*^  have 
devised  techniques  involving  the  measurement  of  the  ratios  of  the  parallel 
and  perpendicular  components  of  the  reflectance.  A  different  technique 
based  on  similar  principles  was  described  recently  by  the  author.*® 

**H.  R.  Philipp  and  G.  A.  Taft,  Phys.  Rev.,  Vol.  113,  p.  1002  (1959). 
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For  a  discussion,  see  F,  Stern  in  Solid  State  Physics,  Vol.  15,  ed. 
by  F.  Seitz  and  D.  Turnbull,  New  York:  Academic  Press  (1963),  p.  300. 

i3R.  J.  Archer,  Phys.  Rev.,  Vol.  110,  p.  354  (1958). 

14 

H.  R.  Philipp  and  H.  Ehrenreich,  Phys,  Rev.,  Vol.  129,  p.  1550 
(1963). 

15 

D.  T.  F.  Marple  and  H,  Ehrenreich,  Phys.  Rev.  Ltrs.,  Vol.  8, 
p.  87  (1962). 

G.  Avery,  Proc.  Phys.  Soc.  (London),  Vol.  B  65,  p.  425  (1952). 

*^R.  E.  Lindquist  and  A.  W.  Ewald,  J.  Opt.  Soc.  Am.,  Vol.  53,  p.  247 
(1963). 

18R.  F.  Potter,  J.  Appl,  Opt.,  Vol.  4,  p.  53  (1965). 


With  this  newer  technique,  a  reflectometer  was  used  to  measure  the 
optical  constants  of  electropolished  germanium  samples  in  the  spectral 
region  from  0.5  to  3.0  eV  at  300*  and  120*K. 

MEASUREMENT  OF  X  *  AND  *n  AND  DETER¬ 
MINATION  OF  THE  OPTICAL  CONSTANTS 

The  experimental  technique,  having  been  fully  described,^  will  not 
be  repeated  here  except  to  say  that  the  spectral  source  is  a  high-intensity 
Bausch  and  Lomb  monochromator  with  a  tungsten  lamp.  Suitable  Corning 
glass  filters  and  interference  filters  are  used  for  separating  orders  at 
the  longer  wavelengths.  The  polarizer  is  a  Glan-Thompson  prism;  the 
detector,  a  lead  sulfide  photoconductor.  The  intrinsic  germanium 
was  kindly  provided  by  T.  Donovan  of  the  Michelson  Laboratory,  Naval 
Ordnance  Test  Station,  China  Lake,  and  had  been  electropolished  by  him. 
The  resolving  power  of  the  monochromator  is  between  60  and  200  over 
the  spectral  range  given  here. 

The  experiment  is  one  of  measuring,  as  functions  of  the  radiation 
frequency,  both  the  ratio  =  R||/Ri  at  the  pseudo-Brewster  angle  (i.e., 
that  angle  for  which  8  is  a  minimum)  and  the  pseudo-Brewster  angle 
(<|>jj)  itself.  The  symbols  R||  and  Redesignate  the  reflected  intensity 
for  the  component  having  the  electric  vector  parallel  and  perpendicular, 
respectively,  to  the  plane  of  incidence;  and  the  pseudo-Brewster  angle 
is  expressed  in  the  quantity  X2  =  sin2  $B  tan2  Figure  2  shows  these 
two  quantities  over  the  spectral  range  at  300 *K,  and  Figure  3  shows  the 
same  quantities  at  120*K.  The  curves  for  Rn  have  been  corrected  by 
subtracting  a  quantity  A(hv)2  where  A  =  4  X  10"^(eV)"2.  This  corresponds 
to  an  oxide  layer,  approximately  20  A  thick  and  having  a  refractive  index 
of  2.0,  on  the  surface  of  the  sample. 

An  examination  of  Figures  2  and  3  demonstrates  the  power  of  this 
type  of  spectroscopy.  One  can  see  certain  spectral  features  present  in 
the  raw  data  prior  to  any  reduction  to  the  optical  quantities.  The  observed 
quantities  have  shifted  both  in  temperature  and  in  value.  The  quantity 
X  2  shows  five  spectral  features  in  this  region,  while  fin  shows  two. 

Their  relationship  to  the  optical  constants  will  be  explained. 

It  was  shown  by  the  author2®  that  the  quantities  X„  and  An  (for  the 
present  range  of  values)  can  be  related  to  each  other  through  the  quan¬ 
tities  a$,  Pjj,  with  an  error  of  less  than  one  percent,  as  expressed  in 
equation  (1). 
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The  quantities  and  €2  for  the  same  two  temperatures  are  shown  in 
Figures  4  and  5.  One  sees  the  marked  degree  of  correspondence  between 
and  the  quantity  €\’,  fin  and  €2  are  not  so  close  in  their  correspond¬ 
ence,  although  similarities  do  exist.  The  quantity  shows  small  but 
quite  definite  spectral  features  in  the  vicinity  of  0.80  and  1.10  eV  at 
300*K,  and  near  0.90  and  1.20  eV  at  120*K.  The  basis  for  this  structure 
is  given  in  Figure  6,  which  shows  the  experimental  determination  of  the 
quantity  Xn^  in  the  energy  region  between  0.7  and  1.30  eV.  The  accuracy 
shown  corresponds  to  a  precision  of  measurement  of  the  pseudo-Brewster 
angle  of  approximately  5  minutes  of  arc  or  less.  The  features  at 
0,80  and  0.90  eV,  as  well  as  those  at  1.10  and  1.20  eV,  are  quite  real 
and  give  rise  to  those  respective  features  in  €}•  In  Figure  7,  the  feature 
at  1.97  eV  at  300*K  and  near  2.03  eV  at  120*K  has  not  been  reported  by 
other  investigators.  The  two  prominent  features  in  €2  data 
2.1  <  hv  <  2.3  eV  correspond  to  the  A  doublet,  which  is  a  prominent 
feature  observed  in  normal  reflectance. 21  »22 


DISCUSSION 

Figures  4  and  5  exhibit  five  spectral  features,  four  of  which  have  been 
reported  by  other  investigators.  In  no  single  experiment  have  all  five 
spectral  features  been  observed.  They  are  all  attributed  to  various  inter 
band  transitions  between  the  valence  band  and  the  conduction  hand  struc¬ 
ture  of  the  germanium  crystal.  The  0.80  eV  transition  is  the  well-known 
direct  gap  at  the  T  point  corresponding  to  the  ^  —  transition, 

21  T_  — 

J.  Tauc,  op.  cit. 

22 

M.  Cardona,  op.  cit. 
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FIGURE  4.  The  Optical  Constants  for  Electropolished  Germanium  at  300  *K, 

Based  on  Data  of  Figure  2 


FIGURE  5.  The  Optical  Constanta  for  Electropolished  Germanium  at  120 

Based  on  Data  of  Figure  3 


Region  0.70-1.30  oV.  Each  feature  is  preceded  by  a  positive  slope. 


ha/  ( ev) 


FIGURE  7.  Detailed  Plots  of  the  Experimental  Quantity  X,^  in  the 
Vicinity  of  2.0  eV.  The  symbols  correspond  to  separate  determina 
tions  and  indicate  the  precision. 


which  has  been  observed  in  transmission  measurements.^  At  120*K  the 
peak  is  approximately  0.38  eV,  a  shift  of  0.08  eV.  The  next  highar  energy 
transition,  appearing  at  approximately  1.095  eV  at  300 'K  and  at  1.22  eV 
at  120*K,  is  attributed  to  the  transitions  from  the  spin-orbit  split  band 

r35'2-  r£i  transition,  which  has  been  observed  by  Hobden.^*  It  would 
appear  from  Figure  6  that  the  direct  edge  has  a  smaller  temperature 
coefficient  than  the  split-off  valence  band  transition. 

The  two  principal  features  in  €3  &t  300*K,at  approximately  2.1  and 

2.3  eV  (Figure  4),  sharpen  considerably  at  120*K  and  shift  to  2.2  and 

2.4  eV  (Figure  5).  These  two  features  have  been  observed  in  all  normal 
reflectance  measurements  since  they  were  first  seen  by  Tauc.  Initially 
they  were  attributed  to  direct  transitions  from  the  L  point  at  the  edge  of 
the  Brillouin  zone,  with  the  valence  band  splitting  corresponding  to  a 
value  between  0.18  and  0.20  eV.  Detailed  band  calculation  by  Brust, 
Philipps,  and  Bassani^  and  by  Brust^  indicated  that  the  joint  density 
of  states  for  transitions  of  this  energy  occur  at  the  A  direction,  which 
is  in  the  <11 1>  direction  but  not  at  the  edge  of  the  zone.  Valence  band 
splitting  would  still  be  of  the  order  of  0.20  eV.  Phillips, ^7  in  his  analy¬ 
sis,  indicates  that  this  transition  (to  be  classed  Mj)  is  associated  with 

a  saddle  point  giving  a  cusp  in  the  joint  density  of  states  at  the  transition 
frequency. 

All  the  transitions  except  the  Mq  transitions  at  the  L  point  have  been 
observed  in  other  experiments^  in  particular  in  the  Franz-Keldysh- 
Seraphin  effect  by  Seraphin  and  coworkers. The  SKS  effect— the 
application  of  an  electric  field  and  the  observation  of  the  induced  changes 
in  the  band  structure  by  reflectance  techniques— was  used  in  an  effort 
to  find  Mq  in  a  search  from  1.8  to  2.0  eV.  It  is  suggested  that  the  effect 
of  the  A— transition  at  2.1  eV  is  so  strong  that  the  SKS  observation 
of  the  Mq  transition  at  1.97  eV  is  obscured.  The  extreme  precision  of 
measuring  Xn^,  it  is  believed,  permits  one  to  observe  the  spectral  feature 
at  1.97  eV  at  300 #K  and  at  2.03  eV  at  120#K  (Figures  7  and  8). 

2^W.  C.  Dash  and  R.  Newman,  Phys,  Rev.,  Vol.  99,  p.  1151  (1955). 

2*M.  V.  Hobden,  J.  Phys,  Chem.  Solids,  Vol.  23,  p.  821  (1962). 

25 

D.  Brust,  J.  C.  Phillips,  and  F.  Bassani,  Phys.  Rev.  Ltrs.,  Vol.  9, 
p.  94  (1962). 

2^D.  Brust,  Phys,  Rev.,  Vol.  134,  p.  A  1337  (1964). 

27J.  C.  Phillips,  Phys.  Rev.,  Vol.  133,  p.  A  452  (1964). 
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B.  O.  Seraphin  and  R.  B.  Hess,  Phys.  Rev,  Ltrs.,  Vol.  14,  p.  138 
(1965). 


79 


The  interband  transitions  occur  at  critical  points  in  the  Brillouin 
zone.  These  critical  points  (Van  Hove  singularities^)  occur  when  the 
relation  -  Ej^i)  =  0  is  satisfied,  as  distinct  changes  in  the 

joint  density  of  states  (J\is)  occur  at  these  points.  Transitions  where 
=  V^Ej^i  =  0  and  the  difference  Ej^  -  E^t  is  a  minimum,  are 
known  as  MQ“tyPe*  Such  conditions  occur  for  Ge  at  the  center  of  the 
zone  (F  point)  and  at  the  <11 1>  edge  (L  point).  A  very  large  joint  density 
of  states  maximum  occurs  between  the  valence  and  conduction  bands  of 
Ge  in  the  <11 1>  direction  (A).  The  energy  band  difference  is  character¬ 
ized  as  a  saddle  function  of  the  type  (see  D.  Brust,  footnote  26). 

The  joint  density  of  states  has  a  maximum  for  the  frequency  Ej^  -  Ej^ji 
at  the  saddle  point. 


Because  the  imaginary  part  of  the  dielectric  constant  is  proportional 
to  the  joint  density  of  states  for  a  transition  (^^J^g)*  it  is  a  very 
important  band  structure  parameter.  Brust  ha3  carried  out  a  density 
of  states  determination  for  different  types  of  Van  Hove  critical  points. 


In  Figure  9  are  illustrated  the  joint  densities  of  states  for  two  of  the 
transition  types.  Figure  9a  represents  the  case  for  minimum  extrema 
in  energy  gap,  i.e.,  at  L  or  T  points  (Mg).  while  Figure  9b  represents 
the  case  for  a  saddle  point  (Mj).  The  contributions  6  to  the  real  and 
imaginary  parts  of  the  dielectric  constant  can  be  shown  to  be  as  follows:-^ 


6e^(u)  =  C  c^fw'Mn  jw'2  -  u2  j  dt*)'  (3a) 

Jtt)  -  A 


6e2(w) 


ft*  +A 

±1  c  «,•!■•>* 

b>*  +  b> 

b)1  -  b) 

J b)  -A 
c 

do' 


(3b) 


where  wc  corresponds  to  the  transition  frequency,  and  wc  +  A  and  uc  -  A 
give  the  frequency  bounds  outside  of  which  the  contributions  6c  i  and  6C2 
ar  :  negligible.  Because  the  logarithms  in  the  integrands  are  large  at 
u>  -*  w',  one  expects  principal  changes  in  €i(w)  and  C2(w)  when  €2  (<*>)  and 
6i(u),  respectively,  have  large  and  abrupt  changes. 


From  Figure  9c  it  is  seen  that  €2(w)(Mq)  Has  a  discontinuity  at  wc. 
Equation  (3a)  indicates  that  in  this  case  6c | (u>)  goes  from  a  ^positive  value 
at  u  <  b)c  to  a  maximum  at  w  wc.  Figure  9d  shows  that  €2  (w)(Mi)  also 
has  a  discontinuity  at  uc  but  it  goes  from  a  positive  value  at  w  <  wc  to  a 


29L.  Van  Hove,  Phys,  Rev.,  Vol.  89,  p.  1189  0953). 
^F.  Stern,  op.  cit.,  p.  332, 
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negative  constant  at  w  >  wc.  Equation  (3a)  indicates  that,  for  negative 
£3,  6€j(w)  has  a  negative  slope  between  a  maximum  value  of  at 
o  &  wc  and  a  minimum  value  at  w  Jt  wc.  Equation  (3b)  provides  the 
corollary  that  a  large  negative  slope  in  corresponds  to  a  positive 
extremum  in  b€2>  A  positive  slope  in  6£j  indicates  that  a  sharp  kink, 
if  not  a  minimum,  occurs  in  ££3  for  the  general  case  Mq  (Figure  9a). 

Although  the  slope  for  the  1.97  eV  peak  (300“K)  in  €j  is  not  clearly 
determined,  it  appears  to  be  due  to  a  weak  transition,  presumably  at 
the  L  point;  hence  the  maximum  in  very  nearly  gi'/ca  the  L3  —  Lj 
transition.  The  detailed  data  for  this  spectral  region  ns-ar  2.0  eV  are 
shown  in  Figures  7  and  8.  On  the  other  hand,  peaks  at  0.80  and  1.09  eV 
appear  to  be  preceded  by  positive  slopes  (Figure  6),  with  maxima  in  6£j 
at  energies  very  nearly  equal  to  those  for  the  assigned  transitions.  The 
negative  slopes  in  £j  corresponding  to  peaks  in  £3  at  2.1  and  2.3  eV  are 
strong  indications  that  direct  transitions  take  place  at  those  frequencies 
of  the  Mj  types,  as  indicated  by  Phillips^  *  and  Brust.^  The  assignments 
for  the  five  spectral  features  are  given  in  Table  1. 


TABLE  1.  Transition  Assignments  (eV) 


Temp- 

Type  MQ 

Type  Mj 

erature 

CK) 

r  3/2-r 

25'  2' 

— r 

25'  2' 

L3"L1 

A3/2-* 

A3  Aj 

A3  A1 

300 

0.80 

1.09 

1.97 

2.11 

2.29 

120 

0.88 

1.24 

2.03 

2.22 

2.41 

It  is  of  interest  to  compare  the  optical  constants  determined  here 
with  some  very  precise  measurements  of  electropolished  germanium 
reflectance  made  by  Donovan,  Ashley,  and  Bennett. ^  Their  measure¬ 
ments  were  made  at  room  temperature  but  with  a  precision  of  better  than 
one-half  of  one-tenth  of  a  percent  in  the  spectral  region  from  1.2  to  4.7 
eV.  The  data  were  tabulated  and  the  results  up  to  3  eV  are  shown  as 
open  circles  in  Figure  10,  along  with  the  solid  curve  showing  the  calcu- 


3lJ.  C.  Phillips,  op.  cit. 

32D.  Brust,  Phys.  Rev.,  Vol.  134,  p.  A  1337  (1964). 

33T.  M.  Donovan,  E.  J.  Ashley,  and  H.  E.  Bennett,  J.  Opt.  Soc.  Am., 
Vol.  53,  p.  1403  (1963). 
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FIGURE  10.  Normal  Reflectance  for  Electropolished  Germanium  at  300 *K  as 
Calculated  From  the  Optical  Constants,  Shown  in  Comparison  With  Precise 
Reflective  Measurements 


lated  normal  reflectance  based  on  the  NOLC  data  of  Figure  4.  Some 
selected  points  based  on  Archer's  ellipsometric  data,^*  which  give  values 
of  the  optical  constants  somewhat  different  from  the  present  data,  are 
also  shown  (closed  circles)  for  comparison.  One  sees  that  the  agree¬ 
ment  is  relatively  good  considering  that  Donovan  et  al.  reported  varia¬ 
tions  of  as  much  as  one-half  of  a  percent  between  different  runs.  It  is 
also  seen  that  the  use  of  normal  reflectance  to  discern  the  Mq  transition 
at  1.98  eV  would  be  extremely  difficult  because  of  the  enormously  steep 
rise  due  to  the  M}  transition  at  2.1  eV.  Normal  reflectance  measure¬ 
ments  of  the  Mq  transition  at  0.8  and  1.1  eV  should  be  possible,  given 
suitable  care  and  precision  such  as  that  prescribed  by  Stern. 

No  absolute  values  of  normal  reflectances  at  low  temperatures  seem 
to  be  available,  although  Cardona  and  Sommers^  show  relative  normal 
reflectance  over  the  spectral  region  for  two  temperatures.  The  appear¬ 
ance  of  the  peaks  at  2.2  and  2.4  eV  (Figure  11)  are  in  good  qualitative 
agreement  with  their  curve  for  90  *K.  It  is  evident  that  the  Mq  peak  at 
2.05  eV  would  still  be  relatively  hard  to  see  using  normal  reflectance. 

As  mentioned  above,  the  normal  reflectance  will  usually  be  representa- 
ti\  e  of  either  the  real  or  imaginary  part  of  the  dielectric  constant  in 
most  regions.  But,  as  seen  here,  what  goes  on  in  the  region  between 
1.5  and  2.0  eV  is  rather  indeterminant.  For  that  reason  the  optical 
constants  directly  determined  with  precision  are  useful  quantities. 


34 

35. 


R.  J.  Archer,  op.  cit. 


M.  Cardona  and  H.  S.  Sommers.  Jr.,  Phys.  Rev.,  Vol.  122,  p.  1382 
(1961). 


FIGURE  IK  Normal  Reflectance  for  Electropolished  Germanium  at  120“K  as 
Calculated  From  the  Optical  Constants  of  Figure  5 


SOLID  STATE  SPECTROSCOPY 


CRYSTAL  FIELD  EFFECTS  FOR  Ce3+  AND 
Yb3f  IN  THE  GARNETS 

by 
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and 

R.  A.  Buchanan 
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The  understanding  of  the  anisotropy  in  the  exchange  interaction  between 
the  rare  earth  and  iron  ions  in  the  garnets  depends  on  detailed  knowledge 
of  the  rare  earth  wave  functions.  Complete  crystal  field  calculations 
have  not  been  possible  in  the  past  because  the  available  data  did  not  suf¬ 
fice  for  a  determination  of  all  nine  crystal  field  parameters  in  addition 
to  the  spin-orbit  interaction.  This  lack  necessitated  the  arbitrary  sup¬ 
pression  of  a  certain  number  of  these  parameters  in  the  attempted  cal¬ 
culations.  The  most  systematic  approach  to  date  was  applied  by 
Hutchings  and  Wolf1  (HW)  to  the  spectrum  of  Yb3+  in  a  number  of  garnets. 
These  workers  used  all  nine  crystal  field  parameters,  imposing  on  them 
a  set  of  constraints  in  which  it  was  required  that  ratios  among  fourth- 
and  among  sixth-order  parameters  conform  to  those  arising  from  a 
detailed  point  charge  computation.  By  adjusting  the  five  remaining 
parameters,  Hutchings  and  Wolf  were  successful  in  obtaining  a  reasona¬ 
ble  fit  for  the  three  J  *  5/2  levels  and  for  two  sets  of  g-factors. 

Since  the  appearance  of  the  above  paper,  the  authors  have  obtained 
additional  experimental  data.  These  data  include  a  new  level  of  Yb3* 
in  yttrium  gallium  garnet  (YGG)  at  —900  cm’1,  a  complete  determination 
of  the  4f  energy  levels  for  the  electronically  analogous  ion  Ce3+  (also  in 
YGG).  and  the  exchange  splittings  of  certain  levels  of  Ce3*  in  YIG.^ 

A  program  has  been  initiated  that  is  designed  to  obtain  as  complete  a 

*M.  T.  Hutchings  and  W.  P.  Wolf.  J.  Chem.  Phys..  Vol.  41.  p.  61? 
(1964). 

*K.  A.  Wickersheim  and  R.  A.  Buchanan,  Fifth  Rare  Earth  Research 
Conference,  Ames,  Iowa,  30  August- 1  September  1965.  The  majority 
of  the  data  mentioned  here  were  available  at  the  time  of  the  Ames  pre¬ 
sentation. 
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set  as  possible  of  crystal  field  and  exchange  parameters  in  the  garnets 
and  to  uncover  any  systematic  patterns  in  the  behavior  of  the  various 
rare  earth  ions. 


As  a  first  step,  the  analysis  of  Ce^+  has  been  undertaken  with  a  view 
towards  testing  the  applicability  of  the  HW  method,  and  an  appropriate 
computer  program  has  been  formulated.  To  facilitate  level  identification 
it  was  found  convenient  to  work  from  the  outset  with  basis  functions 
belonging  to  the  irreducible  representations  of  the  octahedral  group.  In 
this  form  the  non-cubicity  of  the  calculated  wave  functions  can  be  deter¬ 
mined  by  inspection.  As  a  check  on  the  program,  the  spectrum  of  Yb^+ 
in  YGG  was  recalculated  using  the  HW  parameters,  and  found  to  be  in 
agreement  with  their  results.  The  cubic  parentage  of  the  calculated 
wave  functions  was  found  to  be  clear-cut.  The  corresponding  cubic 
labeling,  as  given  in  Figure  1,  follows  a  sequence  Ty,  F^,  Tg  for  the 
J  =  7/2  levels,  rather  than  Tg,  Tg  as  commonly  assumed.  It  is  felt 
that  proper  identification  of  the  quartet  (Tg)  parentage  is  a  generally 
significant  consequence  of  the  re-examination  of  the  HW  calculation,  as 
well  as  being  of  special  importance  in  the  analysis  of  the  analogous  (but 
inverted)  level  scheme  for  Ce^+. 


Next,  a  preliminary  attempt  was  made  to  obtain  a  fit  for  the  Ce^+ 
spectrum,  utilizing  the  HW  parameters  for  Yb^+  modified  by  three 
adjustable  multiplicative  factors  which  are  applied,  respectively,  to  the 
second-,  fourth-,  and  sixth-order  potentials.  Implicit  in  this  approach 
is  the  assumption  that  the  difference  in  crystal  fields  in  the  two  ions 
arises  principally  from  differences  in  radial  wave  functions  and  isotropic 
shielding  factors.  The  success  of  the  procedure  can  be  gauged  from  an 
inspection  of  Figure  2.  The  fit,  while  not  perfect,  is  unmistakable. 

The  parameters  used  (as  defined  by  HW)  are  given  by 


s  -52. 

A*  a  178, 

\°  = 

-347, 

> 

*  0 

=  286, 

*  963. 

a6°  *  101. 

%*■ 

-321, 

A6 

a  1841, 

A,6  *  -326,  and  7/2  X  =  2157  cm 
o 


The  second-order  coefficients  differ  from  those  of  Yb*+  by  a  factor  of 
—0.6.  the  fourth-order  by  —1.8,  and  the  sixth-order  by  —1.4.  The  agree¬ 
ment  between  calculated  and  observed  levels  greatly  increases  confidence 


^Viewed  in  the  framework  of  wave  function  space,  the  calculated 
functions  deviate  2G  degrees  or  less  from  the  cubic  basis  functions. 
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FIGURE  1.  A  Comparison  of  Calculated  and  Observed 
Levels,  in  cm"*,  for  Yb*+  in  Yttrium  Gallium  Garnet. 
Except  for  a  readjustment  of  tha  spin- orbit  parameter 
and  the  labeling  of  cubic  parentages,  the  calculated  results 
duplicate  those  of  Hutchings  and  Wolf.  The  level  near 
900  cm“*  was  observed  in  luminescence. 
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FIGURE  2.  A  Comparison  of  Calculated  and  Observed 
Energy  Levels,  in  cm"*,  for  Ce^  +  in  Yttrium  Gallium 
Garnet,  The  cubic  parameters  used  in  the  cubic-only  fit 
are  not  identical  with  those  used  in  the  full  fit.  The  level 
near  402  cm'*  was  observed  very  recently  in  absorption. 
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in  the  present  identification  of  the  Ce^+  levels.  It  also  reinforces  belief 
in  the  approach  used  by  Hutchings  and  Wolf,  at  least  as  a  good  first 
approximation,  and  at  the  same  time  confirms  their  interpretation  of 
the  optical  spectrum  of  Yb^+  in  the  garnets.^ 

One  should  point  out,  on  the  other  hand,  that  the  energy  levels  are 
sensitive  primarily  to  the  predominant  cubic  terms  in  the  potential,  and 
it  is  therefore  difficult  to  draw  precise  conclusions  concerning  the  non- 
cubic  terms  without  some  knowledge  of  the  eigenfunctions  (e,g.,  from 
g-factors).  In  fact,  the  rather  striking  fit  that  is  obtained  with  a  strictly 
cubic  potential  to  the  Tg,  F7,  and  the  mean  of  the  Tg  components  (see 
Figure  2)  suggests  that  the  immediate  environment  of  Ce^+  is  more 
nearly  cubic  than  that  of  Yb^+.  This  impression  is  reinforced  by  the 
relatively  small  splitting  of  the  Eg  quartets,  and  indirectly  by  the  ex¬ 
change  splitting  observed  in  the  upper  Fy  level  of  Ce^+,  this  splitting 
being  far  more  isotropic  than  for  the  corresponding  (ground)  state  of 
ytterbium. 

Apart  from  the  inversion  of  the  energy  level  scheme,  the  most  con¬ 
spicuous  differences  in  the  spectra  of  Ce^+  and  Yb^+  are  caused  by  the 
fact  that  the  ratio  of  crystal  field  to  spin-orbit  splitting  is  much  larger 
in  the  former  ion.  In  Yb^+  the  7/2  and  5/2  levels  are  well  separated, 
and  for  most  purposes  can  be  regarded  as  uncoupled.  This  is  not  the 
case  for  Ce^+  where  important  contributions  are  made  by  matrix  elements 
connecting  7/2  and  5/2  levels.  There  exists  in  particular  a  very  large 
matrix  element  (predominantly  of  cubic  origin)  which  connects  the  two 
F7  levels,  causing  a  further  shift  (away  from  each  other)  of  more  than 
400  cm'*.  This  contributes  to  the  very  large  separation  between  the 
upper  F7  level  and  the  remaining  7/2  levels  as  compared  with  the  rela¬ 
tively  small  spacings  between  the  5/2  levels. 


The  HW  results  indicate  no  electronic  levels  for  Yb^+  between  the 
ground  state  and  517  cm’*,  whereas  D.  L.  Wood  (J.  Chcm.  Phys.,  Vol. 
39,  p.  1671,  1963)  proposes  two  ytterbium  levels  in  this  range,  and 
J.  A.  Koningstein  proposes  one  (Solid-State  Maser  Research  Report 
[Optical],  Sixth  Quarterly  Report,  dated  14  February  1965,  U.  S.  Army 
Electronic  Command  Contract  Number  DA- 36-039- AMC-02333E). 
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